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ABSTRACT 


The motion and stability of spin stabilized spacecraft with mov- 
able external appendages are treated both analytically and numeri- 
cally. The two basic types of appendages considered are: (1) a 

telescoping type of varying length and (2) a hinged type of fixed 
length whose orientation with respect to the main part of the space- 
craft can vary. Two classes of telescoping appendages are considered: 
(a) .where an end mass is mounted at the end of an (assumed) massless 
boom; and (b) where the appendage is assumed to consist of a uni- 
formly distributed homogeneous mass throughout its length. 

For the telescoping system Eulerian equations of motion are 
developed. During all deployment sequences it is assumed that the 

r 

transverse component of angular momentum is much smaller than the 
component along the major spin axis. Closed form analytical solutions 
•for the time response of the transverse components of angular veloci- 
ties are obtained when the spacecraft hub has a nearly spherical .mass 
distribution. For the more general case, a series solution is = 
obtained and this solution is limited by its radius of convergence. 

The comparison of the different approximate analytical methods with 
numerical integration results are studied. and it is observed that the 
oscillatory nature of the responses of the transverse angular velocity 
components reduces rapidly with faster extension rates. 

As an application for spacecraft rescue and recovery, booms are 
extended along all principal axes to (a) detumble a symmetrical 
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spacecraft, and (b) achieve a desired final spin about one of the 
principal axes. From an application of Lyapunov's second method boom 
extension maneuvers can be determined. Numerical examination of 
detumbling for asymmetrical hubs also is considered. The use of tele- 
scoping systems for detumbling a randomly spinning spacecraft to 
achieve a desired final state in a time optimal manner is studied and 
it is found that simple boom extension maneuvers alone can not be used 

to achieve the desired state in minimum time. 

The equations of motion for the hinged system are developed using 
the Quasi-Langrangian and the general Lagrangian formulation. In this 
formulation there is no restriction on the location of the hinge 
points. 
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time varying coefficients in the approximate 
equations for fu, fi 2 

coefficients in the series solution for h x (t) 

for deployment when ' 3 ' axis is a symmetry 
axis, b(t) = a x (t) « a 2 (t) 

coefficients in the series solution for h 2 (t) 

constants appearing in the approximate analytical 
solutions for h x , h 2 for the case of a nearly 
spherical hub 

boom extension rate 

components of the angular momentum vector along 
the principal axes 

assumed constant value of 113 during nominal 
deployment maneuver 

instantaneous values of principal moments of 
inertia 

hub principal moments of inertia 

moments of inertia at the switching time T3f 
in the recovery sequence to achieve final 
spin about the ' 3 ' axis 

/ Tf dt = If, cost functional for time optimal control 
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K t) 
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m 


= 2pc 3 

= constant length of hinged appendages 
= time varying length of telescopic appendages 
» mass of the main part of the spacecraft 
= boom end mass 
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T 3 f 
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Vi 


V^/cm 

V M /o 


Vnvj/o " n 

V 0 /cm 

X(t) 

itt) 


a 1 ,a 2 


constraint on the control vector such that, 
||u||<m* 

2 me 2 

Constants appearing in the solutions for w i(t) 
for the asymmetrical deployment of booms with 
uniformly distributed mass 

offset of the hinge point(s) from the *3' axis 

ti me 

kinetic energy 

switching time in the recovery sequence to 
achieve final spin about the '3 1 axis 

control vector 

Lyapunov function 

inertial velocity of the i~ mass of the 
(hinged) system 

velocity of the main part of the spacecraft 
with respect to the system center of mass 

velocity of the main part of the spacecraft 
with respect to the center of the coordinate 
system ( 0 ) 

velocity of the mass in the system with 
respect to the center of the coordinate system 

velocity of point ‘o' with respect to the 
system center of mass 

State vector 

solution of controllable norm-invariant 
system under unique time optimal control 
u (t) 

coordinates describing the orientation of the 
hinged appendages relative to the hub 
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1> W 2> 


U). 


ft 

p 
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^ 0 , ^0 


T-j (t) 



= constant value of u> 3 when two pairs of booms are 
extended (symmetrically) parallel to the '3' axis 

= angular velocity components along the principal 
axis 

= desired final value of w 3 (w 3 ^) 

= mass density per unit boom length 

= phase angles appearing in the solutions for 
wi(t)» W 2 (t) and determined from conditions at 
t = 0, t = Irrespectively 

= /b(t)dt 

= elements of control torque vector u 


= indicates time differentiation 
= indicates initial conditions 
= indicates norm of a vector quantity 
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I. INTRODUCTION 


A number of spin stabilized spacecraft have long appendages which 

nominally lie in the plane of rotation perpendicular to the desired 

spin axis. These appendages might be on-board antennas which must be 

extended in orbit after the initial injection sequence. The extension 

of on-board antenna booms is usually done with the use of motors 

located in the central hub of the spacecraft. The dynamical aspects 

1 

of such spacecraft have been discussed in the recent literature. 

Of special interest is the stability of the system during the 
initial extension of boom-type telescoping appendages. An early 
investigation considers this problem for the case of telescoping type 

2 . 

appendages consisting of two end masses at the ends of massless rods. 

It is assumed that the extension maneuver is restricted to a plane 

which is perpendicular to the nominal spin axis and that both the 

system angular momentum and kinetic energy are conserved during this 

maneuver. In addition the transverse components of angular velocities 

are assumed to be zero during extension. Under these assumptions it 

is seen that the resulting Lagrange equations of motion will yield an 

approximate analytical solution for limiting values of initial to 

2 

final moment of inertia ratios. A more recent treatment using an 
Eulerian formulation considers the extension of rigid booms where the 
transverse component of angular momentum remains less than the polar 
component throughout extension. For the special case where the spin 
axis is an axis of symmetry the linearized equations can be solved 
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analytically. A similar approximate solution has also been obtained 

4 

previously under the same type of assumptions. 

The first phase of the current study will examine the three 
dimensional motion of a general spinning spacecraft system with mov- 
able telescoping appendages during the initial deployment manuever. 
During all initial (nominal) deployment sequences, in accordance with 
actual practice, it will be assumed that the transverse component of 
angular momentum is much smaller than the component along the major 
spin axis. The dynamics of such a system will be studied using a 
variety of analytical techniques for special cases and numerical 
methods for the general case. 

It is thought that by using movable and/or extendible booms the 
recovery of a tumbling spacecraft by passive means may be feasible. 
Methods of recovering spinning satellites to a flat-spin condition 
using spin-up thrusters and multiple combinations of thrusters were 
examined in a recent paper. It was concluded that the use of such 
thrusters for the recovery operation are often limited by the weight 
and propellant capacity of the thruster system, and also the reliabi- 
lity problems associated with multiple thrusters in sequence. Kaplan 
describes an alternate recovery system which utilizes a movable-mass 

control device that is internal to the spacecraft and can move along 

6 

a fixed direction track. This device is activated upon initiation 
of tumble and is programmed via a control law to quickly stabilize 
motion about the major principal axis. In a recent related paper it 
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was concluded that the mass track should be placed as far as possible 
from the vehicle center of mass and be oriented parallel to the maxi- 
mum inertia axis; in addition the performance of the control system 
can be improved through larger mass amplitudes along the track and 
also larger mass sizes. 

It is apparent that the location and displacement amplitude of 
any internal control mass will be limited by the physical dimensions 
of the space vehicle. Externally movable appendages could allow for 
a greater range of location and displacement amplitudes of such a 
system; however, as the size of the appendages increases the flexibi- 
lity problems associated with such structures would have to be con- 
sidered. 

Of interest in this study will be the consideration of the 
detumbling dynamics of a spacecraft system with extensible boom- type 
appendages along the principal axes. The recovery maneuver from an 
initial tumble is designed to reach either of two final states. (1) 
close to a zero inertial angular velocity vector and (2) to approxi- 
mate a final spin about a principal axis. (It is thought that small 
terminal residual angular rates could then be removed by temporarily 
activating on-board damping systems). A key advantage of this type 
of system would be its potential reuse for subsequent detumbling 

recovery operations as the need arises. 

For the case of time optimal three-axis control of the nonlinear 
norm invariant system it is an established fact that the control 
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torques about each axis must be proportional to the instantaneous 
angular momentum components about each axis, respectively; it is 
doubtful that such a time-optimal control torque could be generated 
only by boom extension techniques. However it may be possible to 
consider a combination of movable end masses and optimal control jets 
for three-axis control of a tumbling spacecraft. In the event of 
jet failure the movable end masses could certainly be used as a back- 
up re-usable system for detumbling (even if they cannot effect time- 
optimal recovery). 

Other types of spacecraft employ fixed length appendages (hinged 
systems) whose orientation relative to the main spacecraft is changed 
during the deployment maneuver. The dynamics of this type of fixed 
length appendages during the deployment maneuver with rigid appendages 
have been studied only for the case where the transverse components 

of the angular velocity vector are assumed to be zero throughout 

2 . , 
deployment and where the hinge points are located on the hub s pnn- 

2 

cipal transverse axes. The general three dimensional deployment 
dynamics of such a hinged system will be considered here without any 
restriction on the location of the hinge points. 
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II. MOTION DURING DEPLOYMENT OF TELESCOPING SYSTEM . . 

1. General Considerations 

The motion and stability of spin stabilized spacecraft with 
telescoping appendages are studied both analytically and numerically. 
The telescoping appendages considered here are of varying length and 
could represent extensible antennas or a tether connected to the main 
part of spacecraft. Two types of telescoping appendages are consid- 
ered (Fig. 2.1): (a) the case where an end mass is mounted at the 

end of an (assumed) massless boom (end mass moving) and (b) where 
the appendage is assumed to consist of a uniformly distributed, homo- 
geneous mass throughout its length (uniformly distributed mass moving). 

The torque free equations for a spacecraft with varying moments 
of inertia are: 

Hi ~ “ CO^h g 

h 2 - Wj hs — w^hj 

h3 = W2hi - a>ih2 

where h^(t) = I-j ( t) w^( t ) 

Making the approximation: 

IM> IM <<: l h 3 l 

and h 3 = h 0 = const. 


( 2 . 1 ) 

( 2 . 2 ) 


5 



the equations for h* and h 2 become 

• - 

hi = -a 2 (t) h 2 

» 

h 2 = aj{t) hj 

Here a x ( t) and a 2 (t) are defined as: 

(I 3 (t) ' Ii(t)) u 
® l(t) = I 3 (t) hit) h ° 


(2.3) 

(2.4) 


(2.5) 


a2{t) 


(I 3 (t) - I 2 (t)) 
l 3 (t) I 2 (t) 


( 2 . 6 ) 


As a special case when the spin axis is an axis of symmetry (a^t) 
a 2 (t) = b ( t) ) during deployment, Eqs. (2.3) and (2.4) become: 

hi + b ( t ) h2 = 0 (2.7) 

ha - b(t) hj = 0 (2.8) 


Introducing t = /b(t)dt the above equations reduce to, 


dh a 

— + h 2 = 0 (2.9) 

dt 


dh 

- hi = 0 (2.10) 
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The solutions to Eqs. (2.9) and (2.10) can be written as 

h x (t) = q* cos t = q* cos (Aftjdt + #*) (2.11) 

h 0 (t) = q* sin t = q* sin (/ t b(t)dt + $*) (2.12) 

2 o u 0 

The solutions given by Eqs. (2.11) and (2.12) are identical with those 
given in Ref. 3. 

For the general case where there is no axis of symmetry during 
deployment the following approach is taken. Di fferentiating Eq. (2.3) 
with respect to time and using Eq. (2.4), the resulting equation for 
hj is: 

hj - h, + ai(t) a 2 (t) hi = 0 

a 2 (t) 

Similarly the equation for h 2 results as: 
a,(t) * 

h2 ** h2 + ai (t) a2 ( t) h2 “ 0 

Eqs. (2.13) and (2.14) will be extensively used in the following 
sections. 

2. End Mass Moving 

a. Analytical Solution for Asymmetrical Deployment 
The telescoping system, where the end masses are attached at 
the end of massless rods along the ’2* axis, is shown in Fig. 2.1(a). 
The extension of the masses is assumed to originate from the center 


(2.13) 


(2.14) 
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of the spacecraft hub. The moments of inertia about the three prin- 
cipal axes during the extension are 


I, = 


I* + 2ml 2 


I, = I' 


I. = it + 2m l 2 


(2.15) 


3 * 3 

For. a uniform extension rate, l = ct, and introducing P = 2mc 2 Eqs 
(2.15) become: 


* 2 
Ij s I* + Pt 


U = I 


(2.16) 


I, = 


I* + Pt 


From Eq. (2.6) 


2 * 2 * 

llh - I3I2 


a 2 (t) _ 

a 2 (t) I 3 I 2 CI 3 -IZ) 


(2.17) 


Using Eqs. (2.16) and their derivatives, Eq. (2.17) can be written 
as: 


2I*Pt 


a 2 (t) 

a 2 U) = ( I 3 +Pt 2 )( Ig-Ig+Pt^) 


(2.18) 


From Eqs. (2.5) and (2.6) with Eq. (2.17), a^t) a 2 (t) is obtained as 
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(2.19) 


<I 3 - Ii) (Is - I 2 ) 2 
ai(t) a 2 (t) = j 2 j 1 h 0 

312 

o? - it)d* + pt 2 - 1 *) 

= (i* + pt 2 ) 2 ( 1 * + pt 2 ) 1 * 0 


(2.20) 


Eqs. (2.18) and (2.20) are used in Eq. (2.13) to obtain the fol- 
lowing second order ordinary differential equation for hi : 


hi 


2I*Pt hj 

(i* + pt z )a* + 1 * + pt 2 ) + 


df - i?)(i? - n + pt 2 )hgh t s 0 

(i* + pt 2 ) 2 (i* + pt 2 ) 1 * 


( 2 . 21 ) 


From Eq. (2.5), 


^ 1 ( t ) ^ 1 1 1 3 “ 1 3h 

MO = hhih-h) 


( 2 . 22 ) 


The second order differential equation for h 2 is obtained in a procedure 


similar to that used for hi. Using Eqs. (2.22), (2.16) and (2.20) in 


Eq. (2.14) we obtain the result as: 

ii- + (I* + I* + 2Pt 2 ) 2Pt h 2 
(I* + Pt 2 )(I* + Pt 2 ) 

(I* - I*w* - I* + Pt 2 ) h 2 h 2 
+ '{I* + Pt 2 ) 2 (If + Pt 2 ) I* 


(2.23) 
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It can be seen that Eqs. (2.21) and (2.23) contain time dependent 
coefficients with comple'x regular singular points for the case of 
positive boom extension rates (P > 0). 

As a special case, for a nearly spherical hub (I* = ~ I* " I*) 
Eqs. (2.22) and (2.23) reduce to: 


h, - 


21* Pt hj 


!"(!*+ Pt 2 ) (21* + Pt 2 ) 


= 0 


(2.24) 


4Pt h, 


h + 


' 2 (I* + Pt 2 ) 

Eq. (2.24) can be written as: 


- 0 


(2.25) 


in _ 2i *pt 

ha ~ (I* + Pt 2 ) (21* + Pt 2 ) 


* 2Pt 


I* + Pt 2 21* + Pt 2 


(2.26) 


Integrating with respect to time, 


h, - C 


i* + pt 2 
21* + Pt 2 


(2.27) 


Integrating once again, we have, 


hj(t) = C 



t?Tvp 


+ D 


(2.28) 
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where C and D are constants which are determined from the initial 
conditions. From Eqs. (2.2), (2.16) and (2.28), the solution for 
u 1 (t) is given by: 


IDj(t) 


C[t - 



tan' 1 


f2I*/P 


] + D 


I* + Pt' 


(2.29) 


At t = 0, from Eq. (2.29), Uj (0) = D/I* (2.30) 

From Eq. (2.27), h^O) = 1*^(0) = C/2 (2.31) 

Eqs. (2.30) and (2.31) are used in Eq. (2.29) to obtain the final 
expression for ^(t) as: 

ujj(O) + 2w^(0)[t - /l*/2P' tan ^ /2 1 *7 P ^ 

Wj(t) = — " p^2 (2.32) 

1 + 

Eq. (2.25) can be written as. 



4Pt 

(I* + Pt 2 ) 


(2.33) 


This is an exact differential equation and after integrating with 
respect to time, 


h2 


E 

(I* + Pt 2 ) 


(2.34) 
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Integrating once again, we get 


h 2 (t) = pr 


1 


T* I* * I* 3/2 

2~p{t 2 + ~*p) 2(—p) 


' -1, t 
tan ( 


/l*/P 


+ F 


(2.35) 


where E and F are constants which are determined from the initial con- 
ditions. From Eqs. (2.2), (2.16) and (2.35), the solution for o. 2 (t) 

is given by: 


u 2 (t) = p 


2 ^ t 2 + ') 2(^) 3/2 


tan ( /iV 


+ F 


(2.36) 


At t = 0, from Eq. (2.36), w 2 {0) - F/I* (2.37) 

From Eq. (2.34), ( 0) “ I*«2(0) “ E/(I ) (2.38) 

Using Eqs. (2.37) and (2.38) in Eq.(2.36), the expression forw 2 (t) 
becomes 





(2.39) 


In general, for a symmetrical spacecraft, the initial conditions 
^(0) and w 2 (0) can be related to w^O) and u) 2 (o) from the torque-free 
precession before the extension begins. However, it should be noted 
that such initial angular accelerations may also be caused by other 
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types of external perturbations. 

b. Series Solution 

As Eqs. (2,21) and (2.23) cannot be solved in closed form except 
for the special case of a nearly spherical hub, a series solution is 
developed for h,(t) and h 2 (t). A similar type of series solution has 
been previously used to predict the planar librational motion of a 
gravity-gradient satellite during boom deployment. Here t = 0 is 
an ordinary point of Eqs. (2.21) and (2.23), and the radius of con- 
vergence R is the smallest value of* 



The series solution for h 1 may be expanded about t O.in the 


form: 


l *n tn 


(2.40) 


n=o 


Substituting this into Eq. (2.21), we have: 


21 2 Pt l a n *n* t 


n-i 


l a n *n»(n-l )t 
n=o 


n -2 


n=o 


(I* + Pt 2 )(I* + 1* + Pt 2 ) 


+ 


0* - Ipd* - I* + Pt 2 ) h 2 „ 

l ! — i t JL y a t n = o 

(I* + Pt 2 ) 2 (It + Pt 2 ) ^ n = 0 


(2.41 ) 


Rearranging and collecting terms in Eq. (2.41), after a number of 
algebraic manipulations, we get. 


13 



li t— 1 8 


+ 



f [ Ej (n-2)(n-3)-Q,(n-2) +Li]a n _ 2 t n ' 2 
n =2 


n _ 2 

[ Fi (n-4) (n-5) -OjCn-4) + Hi ] a n _ 4 t 


+ l [ Gi (n-6) (n~7) - K x (n-6) + N 3 ] a n _ 6 t n ~ 2 

n-6 


+ H. I (n-8) (n-9) a n _ 8 t"' 2 = 0 

n=8 

where Dj = 1*1* (Ip 2 0 2 + Ip 

E, = 1*P I (Ip 2 Ut + IJ + Ip + 2 Ipf (I* + IP' 

Fl = iV[(ip 2 + ipi: + ip] 

Gj = lV[ it + Ij + 3 1* 1 
1 2 

Hj = i; P 4 ; Q, - 2 I? 1* 3 UP 2 P 
J, - 2dt + Ipdp 2 P 2 s k^poP'p 3 
L x = h 2 (I* - Ip ( (Ip 2 - Op 2 1 

\ = Z K l t (I 3 - ! P P 5 N 1 = h 0 (I t ' T P p2 


(2.42) 
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From Eq. (2.42), 

when n = 0, a Q (0)(-l) = Q ; i.e. a Q f 0 in general; 
when n = 1, D 2 a } (1)(0) = 0 ; I.e. a z f 0, in general; 
.when n = 2, a 2 (2 ) ( 1 ) + L 2 (a Q ) = 0, and 



Si mi l arly. 


(~Qj + L i) a 
Dj-3-2 1 


(2.43) 


(2.44) 


_ _ 1(2*1 - Ej - 2Qj + Lx)a 2 + Hj a 0 ] 


(2.45) 


_ [(3-2-Ej - 3Qx + Li)a 3 + (-J) + Mj) a x ] (2.46) 

a 5 — — — : 

Dj • 5 * 4 

The general form for the first ten terms can be written as (where it 
is understood a ■ « 0 for j < 0): 

■J 

a n = - 'i{n-2)(n-3)E 1 -(n-2)Q 1 + L^} a n _ 2 


+’ {(n-4)(n-5)F 1 -(n-4)J i + M J ) a^ ^ 

+' {(n-6)(n-7)Gi-(n-6)Ki + Nj } a n _ 6 

+ (n-8)(n-9) H 1 a^]/ [Dj -n.(p-l )} (2.47) 
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It is seen that the odd coefficients can be related to aj and the 
even coefficients to a 0 .' The solution for h,(t) is written as 


hj(t) = I a n t' 
n=o 


= a 0 I F j ( t) ] + aj [Gj(t)] 


(2.48) 


where F^t) contains the even coefficients and Gj (t) the odd coeffi- 
cients. The constants a 0 and aj are determined from the initial con- 

ditions as follows: 

a 0 = hj(0) = I* ujCO) 


wj(0) 

(2.49) 

iito) 

(2.50) 


The expression for MjU) is given by: 

hi(t) 


“i 


(t) = 


(2.51) 


I* + Pt 2 


In a similar way, the series solution for ti 2 (t) can be written 


as: 


h 2 (t) = l b n t n 

n=o 


= b [ F (t) ) + b [ G (t) ] (2.52) 

Q 2 1 1 

where F (t) and G ? (t) are similar to Fj(t) and G^t). The constants 
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(2.53) 


are related to the initial conditions according to: 

\ ■ v o) ■ ■; • ! <° | 

- h 2 (0) . 1* : 2 (0) < s - 54 ) 

The solution for u 2 (t) is then given by 


<o 


it) = h jt)/i: 


2 ' ' 2 ' " 2 
The expression for u> 3 (t) results directly from the conservation of 
h 3 C h o ) : 


(2.55) 


co 3 


(t) = 


r + p t 

3 


(2.56) 


c. Numerical Results. 

In this section the results of numerical integration of the non- 
linear differential equations of motion for the most general case are 
presented. The purpose of the numerical investigation is twofold: 
first, to verify the analytical results obtained and, second, to com- 
pare the motion for different cases considered. The numerical 
integration is carried out using the IBM 1130 electronic computer. 

The RKGS and RKSCL subroutines are used to integrate three nonlinear 

12 

equations with time varying coefficients. The subroutine RKGS uses 
the Runge-Kutta method for the solution of initial value problems. 

The purpose of the Runge-Kutta method is to obtain an approximate 
solution of a system of first order ordinary differential equations 
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with given initial values. It is a fourth order integration proce- 
dure which is stable and. self starting; that is, only the functional 
values at a single previous point are required to obtain the func- 
tional values ahead. For this reason it is easy to change the step 
size at any step in the calculations. The entire input of the 
procedure is: (1) lower and upper bound of the integration interval, 

initial increment of the independent variable, upper bound of the 
local truncation error; (2) initial values of the dependent varia- 
bles and weights for the local truncation errors in each component of 
the dependent variables; (3) the number of differential equations in 
the system; (4) as external subroutine sub-programs, the computation 
of the right-hand side of the system of differential equations; for 
flexibility in output, an output subroutine. The subroutine RKSCL 
establishes weighting factors for the error function. 

A typical time response of the components of transverse angular 
velocity for a nearly spherical hub is shown in Figs. 2.2 and 2.3. 

End mass extension rates of c = 4 and c = 1 ft/sec respectively are 
considered where extension is assumed to occur only along the '2' 
axis. For numerical integration Eqs. (2.1) and (2.16) are used to 
obtain the results. The approximate analytical solution given by Eqs. 
(2. 32} and (2.39) and the series solution, given by Eqs. (2.51) and 
(2.55) with ten terms present, are compared with numerical integration 
results. It is observed that the analytical solution corresponds more 
closely with numerical integration results when the extension rate is 
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increased. The series solution can be used only in the initial .part 
of the extension where the analytical solution also gives essentially 
the same result. The series solution is limited by its radius of 
convergence as shown for each case. Fig. 2.4 shows the case of Fig. 
2.3 where the hub is spherical (analytical and numerical integration 
results are the same) and the same initial angular velocity compo- 
nents. Fig. 2.5 is a comparison of the analytical and numerical 
integration results for different initial conditions than those shown 
in Fig. 2.3. It can be seen with the numerical integration results 
that even though the final magnitudes of the angular velocities are 
small, the responses differ predominantly for the intermediate time 
ranges. 

In Figs. 2.6(a) and (b), the effect of varying the hub spin- 
axis moment of inertia is shown using numerical integration. It is 
observed that when the hub spin axis moment of inertia (I*) is 

increased from a spherical one (I 3 = I* = I 2 = 5 slug-ft ), the trans 
verse angular responses tend to become more oscillatory in nature 
during the full deployment time. 

The effect of varying the end mass is considered next and it is 
seen (Fig. 2.7) that the transverse angular velocity amplitudes tend 
to decay more rapidly when the end mass is increased. This type of 
response is due to the increase in moment of inertia when the end 
mass is increased. 
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3. Uniformly Distributed Mass Moving 


a. Analytical Solution for Asymmetrical Deployment 
The telescoping system with the uniformly distributed mass mov- 
ing along the '2' axis is shown in Fig. 2.1(b). The extension of the 
telescoping system is considered to originate from the center of the 
spacecraft hub. The moments of inertia during this deployment can be 
expressed by, 


1 1 = i r + 




i, = r 


(2.57) 


* o ? 
I = T + £ or 
3 3 3 


where P = linear density = mass/unit length. 

With t = ct and K = 2 pc 3 > Eq. (2.57) can be written as 


h = + | t 3 


r I* (2.58) 

I 3 = I* + | t 3 

Following the same procedure as used for the case of the moving end 
mass the angular momentum equations for hj and h 2 , from Eqs. (2.13) and 
(2.14), can be developed to yield the following: 
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(2.59) 


I* Kt 2 


Kt' 


dr I*)d3 - 12 +T-) 


hi- 


,„*.Kt'w,* ,*,W\ ■ ( T*| -1 f T * . M. .Jl* 

(l3+-j-)(I 3 -I 2 +^-) (I 3 + — ) Ui + 3 Jl 2 


.* Kt 

.T 

.3 


h,+ 


^ r T* , Kt 3 \ 2/ T* . Kt 


3 h O h l 


= 0 


h 2 + 


0; ♦ 1; * ^T)Kt 2 . a* - i*hi*-i*.+ t-) 2 

Kt h = h? 


Kt 


Kt 


V 


Kt 3 . 2 . * 


» 0 


. x r\. u , . x f-r . - \ r 

(I 3 + - 3“)(Ii+ 3 ) ( J 3 + 3 ) dl + 3 ) J 2 


(2.60) 


As a special case, for the nearly spherical hub 
(I* =!*=!*= I*), the above equations can be approximated by: 


3 I*hj 


(I* + Ml)t 
3 


= 0 


(2.61) 


2Kt z h 


hz + (i* + Mi) 

3 


= 0 


(2.62) 


Eq. (2.61) can be written as 


hi 

r 


3a 
T77T 


where 


(a 3 + t*)t 


31 


(2.63) 


(2.64) 


Integrating Eq. (2.63) 


hi - R 


t 3 

1 a 3 + t 3 


(2.65) 


where R x is a constant. Integrating again and introducing the initial 
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conditions the solution for wi(t) can be obtained as: 


Ri a . . t a ; *- a, -i 
U1 (0) +T l[t- -zm^^y^un (j^) + ? )] 

_____ 


(t + a) 2 a r -i,2t-a 


Wj(t) = 


( 2 . 66 ) 


The constant Rj cannot be determined from the initial conditions using 
Eq. (2.66), since at t = 0, Ri is indeterminate. A series solution 
(about the ordinary point t - 0) of Eq. (2.61) can be developed to 
yield: 


<Vt) 


wi(0) + R 2 t 4 [ 1 - 
1 + t 3 /a 


4 t 3 
7 P" 


+ 


_4 t 6 
10 F 5 " 


] 


(2.67) 


Here also the constant R 2 cannot be determined from the initial condi- 
tions using Eq. (2.67). 

A different, approach is now adopted using the parent equations of 

Eqs. (2.59), (2.60), i.e. Eqs. (2.3) and (2.4), to evaluate R } in Eq. 

* 

(2,66). The equation for h x is obtained from Eqs. (2.3), (2.6), and 
(2.58) as below: 



( 2 . 68 ) 


For the case of a nearly spherical hub, Eq. (2.68) reduces to the 
following equation, using Eq. (2.64), 
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(2.70) 


Equating Eqs. (2.65) and (2.69) we obtain: 

R r - h o h 2 z 1 * 

Using the initial condition, h2(0), and recalling that h Q = I* w 3 (0) = 
(0) f the constant may be evaluated by, 

Rj = - I* « 3 (0) «, 2 (0) ( 2 - 71 ) 

Using Eq. (2.71) in Eq. (2.66), the solution, Wj(t), results as: 


b) 


0) 


i(t)= 


, (0)-u 3 (0)u,(0) [t- |£n { iLiliT-} - 

A O n + 2 4- _L ^ 2 


t z -at+a 2 /3 


_A{tan 1 ( ^ ^ p )+l}3 


1 + t 3 /a 3 


Eq. (2.62) can be written as 


(2.72) 


^2 = - 
h2 


2 Kt 2 


(i* + ^-) 


(2.73) 


Integrating Eq. (2.73), we get 


bs = — (2-74) 

(I* + Kti) 2 

where Si is a constant. 

Integrating again, and introducing the initial conditions, the solution 
for w 2 (t) results as: 



« 2 (t) = w 2 (0} + 


(a + t) 2 
a 2 - at + t 2 


} 


^2(0) 

3 


a 3 t 

a 3 + t 3 



In { 



,, -l ,2t 
{tan ( 


/r a 


) + 



(2.75) 


b. Numerical Results 

Figs. (2.8) and (2.9) represent a typical comparison of analyti- 
cal with numerical results. Extension rates c = 4 and 1 ft/sec, 
respectively, are assumed for an asymmetrical- deployment only along 
the '2' axis. The analytical result is obtained using Eqs. (2.72) and 
(2.75). Numerical integration results are obtained using Eqs. (2.1) 
and (2.58). The analytical approximation improves with the faster 
deployment rate for the case of the nearly spherical hub. (The same 
type of improvement with faster deployment has previously been noted 
for the case of the moving end mass, Fig. 2.2). Because of the very 
limited radius of convergence of the series solution, a comparison 
with this method of solution was not performed for the case of uni- 
formly distributed mass along the boom. 

The response of both types of telescoping systems, when the uni- 
formly distributed mass is replaced by an equivalent end mass, is 
shown in Figs. 2.10(a) and (b). It is seen that initially both types 
yield approximately the same responses but the amplitudes of trans- 
verse angular velocities are rapidly reduced for the end mass system. 
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This shows the effect of increased moment of inertia, as time 
increases, due to the total mass being placed at the end of the boom. 

The total computer time for numerical integration of the general 
torque free equations with step size of At = 1 sec varied from 135 to 
145 secs for both the moving end mass and uniformly distributed mass 
cases. The extension rates considered were c = 4 and 1 ft/sec for a 
total boom length of 60 ft. The computer time for the evaluation of 
each analytical solution for the above cases considered was about 20 
to 25 secs. 
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■ a. END MASS MOVING 



ORIGINAL PAGE IS 
OP POOR QUALITY* 

b. UNIFORMLY DISTRIBUTED MASS MOVING 


FIG. 2.1. TWO TYPES OF TELESCOPING APPENDAGES 
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FIG. 2.2. COMPARISON OF ANALYTICAL, NUMERICAL INTEGRATION AND SERIES SOLUTION RESULTS FOR A 
NEARLY SPHERICAL HUB WITH EXTENSION RATE c = 4 ft/sec, (END MASS MOVING) 
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FIG. 2.4. COMPARISON OF ANALYTICAL AND NUMERICAL INTEGRATION RESULTS FOR A SPHERICAL HUB WITH EXTENSION 
RATE c = 1 ft/sec. (END MASS MOVING) 



EXTENSION RATE c= 1 ft/sec WITH THE- INITIAL CONDITIONS DIFFERENT FROM FIG. 2.3. (END MASS MOVING) 
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FIG. 2.6(a). EFFECT OF VARIATION OF 1$ FOR FIXED INITIAL CONDITIONS WITH EXTENSION RATE c = 1 ft/sec 
(END MASS MOVING)-NUMERICAL INTEGRATION. 
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FIG. 2.6(b). EFFECT OF VARIATION OF If FOR FIXED INITIAL CONDITIONS WITH EXTENSION RATE c-». 1 ft/sec 
(END MASS MOVING) - NUMERICAL INTEGRATION. 
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FIG. 2.9. COMPARISON OF ANALYTICAL AND NUMERICAL INTEGRATION RESULTS WITH 
EXTENSION RATE c = 1 ft/sec (UNIFORMLY DISTRIBUTED MASS MOVING). 
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FIG 2.10(b) . COMPARISON OF BOTH TYPES Op TELESCOPING SYSTEMS WITH 
EXTENSION RATE c - 1 ft/sec - NUMERICAL INTEGRATION. 



III. USE OF TELESCOPING SYSTEM FOR DETUMBLING 


1 . General Considerations 

The dynamics of detumbling a randomly spinning spacecraft using 
externally mounted, movable telescoping appendages are studied both 
analytically and numerically. The appendages considered are of vary- 
ing length and could represent extensible booms or a tether connected 
to the main part of the spacecraft. Two types of telescoping append- 
ages are considered: (a) the case where an end mass is mounted at 

the end of an assumed massless member (end mass moving) as shown in 
Fig. 3.1; and (b) where the appendage is assumed to consist of a uni- 
formly distributed, homogeneous mass throughout its length {uniformly 
distributed mass moving). 

The extensible boom type appendages are assumed to originate from 
the center of the hub along the three principal axes. The desired 
final states of the system considered are: (1) zero inertial angular 
velocity vector and (2) a final spin about one of the principal axes. 
The necessary conditions for asymptotic stability during the detum- 
bling sequences are determined using Lyapunov's second method. 

2. End Mass Moving 

a. Development of Kinetic Energy 

The configuration of the system, where the end masses are assumed 
to be attached to the end of massless rods along all three principal 
axes is shown in Fig. 3.1. The end masses are assumed to be identical 
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(i.e. m i = m) . The rotational kinetic energy of the system can be 
developed as: 

T = I [{I* + 2m Ul + -£ 3 )) + {I 2 + “2 

+' {I* + 2m(£j + l\)) ..3 + 2m(£j + l\ + l\)) ( 3 -D 

Defining, Ij = I* + 2m {l\ + l\) 

•I 2 = I| + Zm(t\ + l\) (3 * 2) 

I 3 = I* + 2m (zl + l\) , 

Eq. (3.1) can be rewritten as: 

T = 1 [ IjOJj + I 2 o >2 + I 3 W 3 + 2m (l\ + l 2 + l\)\ (3-3) 

If the extension rates are assumed to be constant) Eq. (3.3) can 
be expressed: 

T = l [ I,W* + I,<4 + I,u>*) + non-negative const. (3,4) 

Here the moments of inertia are time varying as the length of the 
booms varies during extension. 

b. Achieve Zero Inertial Angular Rate 
1. Lyapunov Function-Kinetic Energy 

The desired final state of the system is w-j = 0. A suitable 
Lyapunov function, in the state variables w 2 and is the system 
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rotational kinetic energy which can be written as: 

V = T = [ IjWj + I 2 u > 2 + 1 + non-negative const. (3.5) 

The Lyapunov function, V, is positive definite in the state variables 
selected^ for asymptotic stability V will now be examined. 

Differentiating Eq. (3.5) with respect to time, there results: 

V = + I 2 u >2 + + 2 1 1 co 1 oj x + 2I 2 o) 2 cu 2 + 2I 3 u> 3 id 3 ) (3.6) 

The equations of motion can be written, from Eq. (2.1), in the follow- 
ing form: 


• 

hj = l » 3 h 2 - co 2 h 3 = Ijoij + I j»j 

(3.7a) 

* • ■ * 

h 2 * “l h 3 ' “3 h l = 1 2*2 + l 2*2 

(3.7b) 

« • • 

^ 3 = 63 2^'l ” W 1^2 ^ 3 tJ 3 + ^ 3^3 

(3.7c) 

Multiplying Eq. (3.7a) by w 19 Eq. (3.7b) by &> 2 , and Eq. 

(3.7c) by w 3 , 


•we obtain the following: 

((0 3 h 2 - a) 2 h 3 ) Uj 

(iOjhj - Ujhj ) w 2 

(Ug^i “ a} 1^2 ) ^3 
Eqs. (3.8) can be combined to yield: 




1 2*2*2 




; 2 

Vi 


* 2 

V 2 


; 2 

^ 3 W 3 


(3.8) 


Vl“l + ^ 2 * 2*2 + ‘ 


- (IjUj + i 2 (0 2 + I3M3) 


(3.9) 
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Substituting Eq. (3.9) into Eq. (3.6), we obtain 


v = - - (ijWj + i 


A + 1 


A*) 


(3.10) 


From (3.10), we conclude that V is a negative definite function in the 

* * * 

state variables only if 1 2 , I 2 , I 3 > 0. 

Here it is seen that when the rotational kinetic energy is used 
as a Lyapunov function expressed in terms of the inertial angular 
velocity components, that the necessary conditions for asymptotic 
stability are satisfied for positive constant boom extension rates and 
three orthogonally mounted sets of booms along the hub principal axes. 
This means that as time becomes extremely large (and boom lengths 
become infinite) it would be theoretically possible to de-spin a tum- 
bling spacecraft and achieve a zero inertial angular velocity state. 

(Of course, such a situation will, in practice, not occur due to finite 


length appendages, but it will be of interest to simulate how much of a 
random tumble could be removed by this process. The selection of rota- 
tional kinetic energy as a Lyapunov function has also been used by 

7 

Edwards and Kaplan for the system treated in Ref. 6.) 


2. Analytical Solution 

The solutions for the angular momentum of a symmetrical spacecraft 
(Ii « I 2 5 I) during deployment are obtained from Section II. 1 as: 

h,(t) = q* cos (/* b(t)dt + **) (3.11) 

1 0 0 U 
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h 2 U) = q* sin (/* b(t)dt + 4>*) 


(3.12) 



h 3 (t) = h Q = const 

(3.13) 

where 

b(t) = h 

I(t) I 3 (t) 0 

(3.14) 


The moments of inertia about the principal axes are given by 
(Fig. 3.1): 


Ij = I* + 4ml 2 = 

I* + 2Pt 2 . 


1 2 = I* + 4ml 2 = 

I* + 2Pt 2 

(3.15) 

I 3 = I* + 4ml 2 - 

I* + 2 Pt 


Using Eq. (3.15) in Eq. (3.14) 

we obtain: 



h 0 . 1 1 

” 2P { d !“T”F “ d^ + t 2 


} 


(3.16) 


where dj = /I*/2P and d 2 - /I 3 /2P (3.17) 


Introducing Eq. (3.16) in Eqs. (3.11), (3.12), and (3.13) and after 
performing the integration, the solutionsfor the angular velocities 
are obtained as: 


«i(t) 


* ho 1 _ i t 1 , t * 

q 0 cost 2 p { tan dj tan ^ }+ * 0 ] 


I* + 2Pt 2 


(3.18) 
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(3.19) 


.* , «*I- r _L tan" 1 -A - 


<\o sin [ 2 P { 


1 it * 

— tan" — } + ] 

1 2 _J_2 




(t) - 


I* + 2 Pt 2 


I 3 -3(°) 

I* + 2Pt 2 


(3.20) 


where q* and are determined from the initial conditions. Vie observe 
here for large values of t, the solutions for the angular velocities 
lead to the form: 


w.(t) = const/(I* + 2Pt 2 ) , i s 


(3.21) 


This equation indicates that the magnitudes of the angular velocities 
decrease during extension of the appendages, with the square of the 
elapsed time. 


3. Numerical Results 

A typical detumbling maneuver for an initially slowly tumbling 
spacecraft is illustrated in Figs. 3.2 and 3.3. In this example 
because of symmetry the uncontrolled torque-free motion (Fig. 3.2) can 
be theoretically predicted. With an extension rate of 4 ft/sec after 
60 ft, of extension along all three principal axes the angular veloc- 
ity components have been reduced by more than a factor of 10 and, if 
240 ft. of boom could be extended, by a factor of over 300, to a value 
comparable with the orbital angular rate. Removal of this residual 
angular velocity could then be achieved by activating on-board 
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damping devices. 

When the initial angular velocity is increased by an order of 
magnitude the uncontrolled situation (Pig. 3.4) can be recovered as 
shown in Fig. 3.5. For the same extension rate and end-masses the 
order of magnitude reduction in the total angular velocity vector is 
similar to that shown in the slow tumbling case. 

The effect of extension rate on detumbling is illustrated in 
Figs. 3.6(a) - (c). For small extension rates (up to 1 ft/sec) the 
oscillatory nature of the transverse motion is not removed until 
after the first cycle; the advantage of considering higher extension 
rates (at the expense of on-board power) for an initial fast tumbling 
is apparent. It should be noted that at a given time in these figures 
different boom lengths are represented according to the extension rate. 

Numerical examination of other cases for asymmetrical hubs also 
verifies the practicality of using movable appendages for the initial 
detumbling of randomly spinning spacecraft. (Figs. 3.7(a) and (b)). 

The numerical simulation results for an asymmetrical spacecraft are 
compared with the closed form solution for a symmetrical extension and 
it is observed that the closed form solutions are only applicable when 
the asymmetry is small. 

c. Achieve Final Spin About One of the Principal Axes 
1. Lyapunov Function-Modified Kinetic Energy 

The desired final state of the system is; = 0, w 2 = 0 and 
w 3 ~ w 3 f = Using the state variables W|,.w 2 , and u 3 - u> 
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the Lyapunov function is defined as the modified rotational kinetic 
energy, which can be written as: 


V = 1 [ Ijuf + I 2 co| + I 3 (u 3 - a) 2 ] 

Here V is positive definite in the state variables selected, 
entiating Eq. (3.22) with respect to time, we get 

v = 1[ h“i + ' l 2 w l + *3 (“3 " n ) 2 


+ 2IjO)lwi + 2 I 2^2 w 2 + 2 I 3 ( 0)3 - fi ) w 3 ] 

Using Eq. (3.9) in Eq. (3.23), we obtain, 

V = + ^ 2 W 2 + ^3^3^ + 2^ 2l 3 " ^3 W 3 * ^ 3^3 ^ 


For symmetry about the '3' axis during extension: 

• • 

h 3 = ^“3 + 1 3 W 3 = 0 

Eq. (3.25) is used in Eq. (3.24) to obtain: 

V = - \ [ Vl + *2 U 1 + ^ 3 (w 3 ‘ ^ ] 

After rewriting Eq. (3.26) in terms of the state variables, 

V = - 1- [ IjOlj + ^ 2 U 2 + ^ 3 ( u 3 " 1 3^1 (“3 " ii) 

Also from Eq. (3.25), the solution for u 3 (t) is given by, 

Mj(t) = I 3 l0 3 ( 0)/ I 3 


(3.22) 
Differ- 

(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 
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We conclude from Eq. (3.27) that V is negative definite in the -state 
variabl es 

♦ » • 

only if w3 * ft, Ii, I 2 > °> and ! 3 > 0 for w 3 > o 

Thus for the case where a spin about one of the principal axes 
is a desired final condition, a modified form of the kinetic energy 
can be used as a Lyapunov function. Here the final state can be 
achieved by extending all telescoping booms until the desired spin 
rate is reached and then continuing the extension of the set of booms 
along the nominal spin axis until the transverse components of angu- 
lar velocity reach an acceptably small amplitude (within the limita- 
tions of boom length). It should be noted that if we allow w 3 < ft and 
I 3 / 0, there will be a difference in sign between the third and 
fourth terms in Eq. (3.27). 

2. Analytical Solution 

The time at which u> 3 = = ft will be. denoted by T 3 ^. 

At t = T 3 ^, 


I* + 2P (T 3f ) 2 

(3.29) 

it + 2P (T 3f ) 2 

(3.30) 


For t < T 3 ^, the solutions for the angular velocities can be obtained 
from Eqs. (3.18), (3.19) and (3.20). 
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For t > T 3 f , 

(3.31) 

(3.32) 

where If = I* + P(T 3 f) 2 (3.33) 

From Eq. (3.14), and using Eqs. (3.31) and (3.33), we obtain 


= const = I 3 + 2P (t 3 ^) 

I* = I* + P(T 3 f ) 2+ Pt 2 = If + Pt 2 


b(t) = 


(.0 


3 f 


{ 


-* 

Id -1} 


(3.34) 


If + Pt 2 


Introducing Eq. (3.34) into Eqs. (3.11) and (3.12), the solutions for 
the angular velocities for t > T 3 ^ are, 

T 

£f 
T 
f' 


"* f ‘3f 

%«s[ U3f [7p^- {tan' 1 ( 7 p^-)-tan' 1 (7^h-t+T 3f ]+'i' 0 ] 


W l(t) = 


( I jT + Pt 2 ) 


(3.35) 


-* 

q 0 Sin[ “3f [ A 




t 3 f 

(tan-i ( ==) -tan“ 1 ( — 5 — ) }-t+T 3 ] + « 0 ] 

SlpP IpP f 


(I* + Pt 2 ) 


(3.36) 


and from Eq. (3.13), 


w 3 (t) = = const 


(3.37) 


47 



Here q 0 and i|> 0 are to be determined from Eqs. (3.18) and (3.19) at 
t - T 3 ^ and should not be confused with q* and which are determined 
at t s 0. 

For large values of t, Fqs. (3.35) and (3.36) reduce to the form, 

q n cos (u,, x const x t + const) 
u , = — (3.38) 


W 2 


qo sin (w^ x const x t + const) 
I* + Pt 2 


(3.39) 


The above two equations indicate that the frequency of oscillation 
approaches a constant value and the magnitude of the oscillation 
decreases with the square of the elapsed time. 

The time t = T 3f at which the extension of the booms along T 
and '2' axes are stopped can be determined from h3 = 0, yielding the 
result: 


T3f 


_1 

2c 



W3(0) - to3 f 

( 1 ) 


to3. 


(3.40) 


3. Numerical Results 

Figs. 3.8 and 3.9, with extension rates c = 4 and c = 1 ft/sec, 
respectively, illustrate a recovery maneuver which would result in a 
final spin about the *3' body axis with a small transverse residual. 
The booms are extended so that the modified rotational kinetic energy 
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is positive definite and its total time derivative is negative defi- 
nite during the maneuver. All booms are extended until T 3 .p at which 
time. o ) 3 = o) 3 f . Then, only booms along the ± '3' axis are extended to 

reduce the transverse residual components. 

A comparison of the recovery maneuver of an asymmetrical space- 
craft with that of a symmetrical spacecraft to achieve a final spin 
along the * 3 r axis is shown in Figs. 3.10(a) and (b). The calculated 
T 3 f for the symmetrical spacecraft is used for stopping the booms 
along the 'V and 1 2 1 axes. It is observed that using this logic the 
final reaches a lower value ( 1.8 rad/sec) when compared with the 

df 

desired final value (2.0 rad/sec). Also we notice from Fig. 3.10(a), 
the response of ^(t) for the asymmetrical case differs from that of 
the symmetrical case. This is due to the increase in the order of 
the system equations for the asymmetrical extension (i.e. - three 
first order differential equations must now be considered). It should 
be pointed out that after T 3f , for the asymmetrical case, the time 
response of 0)3 is not exactly a straight line as apparently indicated 
in Fig. 3.10(a) but also consists of small amplitude oscillations 
superimposed about this straight line solution. For larger asymmet- 
ries this oscillation would become apparent within the plotting scale 
shown and the difference between w 3 ^ achieved and desired would also 
increase using the open loop control logic of switching the extension 
sequence at a pre-set T 3 ^. 
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3. Uniformly Distributed Mass Moving 


a. Achieve Zero Inertial Angular Rate 
1. Analytical Solution 

The desired final state of the system is ^ =0. The booms con- 
sidered are assumed to have a uniformly distributed mass along their 
lengths. The same procedure as adopted in the case of end mass moving 
can be applied here to obtain the solutions for the angular velocities. 
Here we present only the final results. 

The solutions for the angular velocities are given by: 


q* cos { / t b(t) dt + ip* } 


"i(t)‘ = 


I* + | Kt 3 


(3.41) 




q* sin { / t b(t) dt + } 

o o 


I* + 4 - Kt 3 


(3.42) 


“ 3 ^ = 


I* “3(0) 
I* + | Kt 3 


(3.43) 




3 “ 3 

2 


df a 


2t - d 3 1 (d 4 + t) 

tan "' { “ 6d? lU f d?-d,.t + t 2 } " d 2 /3 


‘4 “4 


1 . ,2t - d4 

— tan -1 { — 

d. / 3 
4 


}] 


(3.44) 
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(3.45) 


d 3 = and d 3 = 21— 

3 2K 4 2K 


b. Achieve Final Spin About One of the Principal Axes 
1. Analytical Solution 

The desired final state of the system is toi = 0, ^2 = 0 an ^ 
w 3 = » 3f = ft. For t < T 3f , the solutions for the angular velocities 

can be obtained from Eqs. (3.41), (3.42) and (3.43). For t > T 3 ^, 
the solutions for the angular velocities can be obtained as: 


Wj(t) 


q 0 cos{ /*■ b(t) dt + ) 


3 f 


I* - + 1 K t 3 
f 3 


(3.46) 


q 0 si h { b(t) dt + ) 

T 3f 

w 2 (t) = 7 ~ (3.47) 

L + 4 K t 5 


“ 3 (0 = “3 f = 

const 

(3.48) 

where: I* = I* + $ 

f ° 

(T 3f ) 3 

(3.49) 


f l b(t)dt = w 
T 3f 


3 I? 1 


K 


6d| 

0 


In { 


(ds + t) 4? 
dl -d 5 t + t' 
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(3.50) 


+ — 1 — tarW f 2t - tis -t + «T i 

d£ A 3f 


i*= i; + |k(t 3 ) 


d»-' 3Jf 

-> 1 / 


(3.51) 

(3.52) 


The time T 3 ^, at which the booms along the 1 1 1 and '2' axes are stopped, 
can be obtained as: 


t = r 3 13 , w 3(°) - w 3f A ! 

3 f ' l 2K 5 


1/3 


(3.53) 
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FIG 31. SYSTEM GEOMETRY FOR END MASS EXTENSION MANEUVER USED 
TO RECOVER A TUMBLING SPACECRAFT. 
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ONTROLLED MOTION - SLOW TUMBLING. 
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FIG. 3.5. RECOVERY DYNAMICS FOR INITIAL FAST TUMBLING 
WITH EXTENSION RATE c, = 4 ft/sec. 
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COMPARISON OF RECOVERY DYNAMICS OF ASYMMETRICAL SPACECRAFT 
WITH SYMMETRICAL SPACECRAFT (EXTENSION RATE c.j * 4 ft/sec). 
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FIG. 3.8. RECOVERY MANEUVER TO ACHIEVE FINAL SPIN ALONG '3' AXIS 
WITH EXTENSION RATE c i = 4 ft/sec 




FIG. 3.9. RECOVERY MANEUVER TO ACHIEVE FINAL SPIN ALONG 
1 3* AXIS WITH EXTENSION RATE q = 1 ft/sec. 
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FIG. 3.10(b). COMPARISON OF RECOVERY MANEUVER OF ASYMMETRICAL SPACECRAFT WITH SYMMETRICAL 
SPACECRAFT TO ACHIEVE FINAL SPIN ALONG '3' AXIS (EXTENSION RATE ci=4 ft/sec). 


IV. TIME OPTIMAL CONTROL 


1 . Combination of Booms and Control Jets 
• (Norm Invariant Principle) 

In this section, we shall consider the control of a norm- 
invariant system which has the property that the Euclidean length of 
the state vector is constant when the control is zero. Here we state 
the problem and the control law to achieve the time optimal control 
of the system from Ref. 9. 

Problem : 

Given the controllable norm-invariant system 

X(t) = g [ X(t); t ] + u ( t ) ; X(0) - 1 (4.1) 

Assume that the dimension of u(t) is the same as the dimension of X(t) 
and that || u(t) | | < m* for all t. Then determine the control which 

forces the system (4. 1 )from the initial state 1 to 0 and which mini- 
mizes the cost functional 

J = dt = T-f, Tf - free (4.2) 

D 

Control Law : 

The unique time optimal control u*(t) that is, the control which 
minimizes the cost J of Eq. (4.2) is given by, 

X*( t ) 

U*(t) = -m* -jj s * (t) || (4-3) 

where X*(t) is the solution of Eq. (4.1) with u(t) = u*(t). The 
minimum value of 0* of the cost J, that is, the minimum time t*. 
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required to force l to 0 is given 

j* = t * = 1 1 « 1 1 (4. 4) 

m* 

The above theory, which deals with the time optimal control of 
norm-invariant systems, is applied to the case of an unsymmetrical 
spacecraft tumbling in torque-free space. The angular momentum 
equations given by Eqs.(2.1) can be rewritten as: 




h?.h3 



( T7‘ TT* hA 


(4.5) 


h 3 e ~ T” ) h i h 2 
i 2 

• 

we find that, -Jl||h(t)|| = M 5.lL = 0 (4.6) 

dt | | h( t) | | 

and so the system represented by Eq. (2.1) is norm-invariant. 

Eqs. (2.1) describe the motion of the spacecraft in the absence 
of any applied torques. A torque vector, u(t), can be generated by 
means of gas jets, reaction wheels, gravity-gradient arrangements, 
etc. At any rate, if u(t) is a control torque, whose components are 
T.j(t), i = 1,2,3, the equations of motion become: 
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h l (t) = [ I 3 (t) ‘ I 2 Ct) ] h 2 t ^ 


L(t) = [ 1 L_] h 3 (t) hj(t) + T 2 (t) (4.7) 

Ii(t) I 3 (t) 

MO = l TKiy- iJtT 1 h i (t) h 2 (t) + T3(t) 

We can immediately conclude that, if the constraints on the control 


torque u(t) are of the form s 


||u(t)|| = /T 3 (t) + T 2(t) + TjltT s m* 

the torque components for time optimal control are: 


m*Ii(t)u>i(t) 


’■ (tl = - -iiiiili 


(4.8) 


T2 (t) = - , * 1 2 (t)l, 2 (t) 
I |h(t) | | 


(t) = - 


m*I 3 (t)u 3 (t) 

||h(t)|| 


(4.9) 


where 1 1 h ( t ) is defined as: 

||h(t)|| = /lj(t) lOj(t) + l 2 (t) W 2 (t) + I 3 (t) u> 3 (t) (4.10) 

This means that, in order to reduce the angular momentum vector R(t) 
to zero in the shortest possible time, the torque vector u(t) must 
point in the opposite direction to the angular momentum vector h(t) 
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9 

and the torque U(t) must be as large as possible. 

For the case of symmetrical spacecraft ( I i ( t ) = l 2 (0)s the tor- 
que components required to reduce the transverse angular momentum to 
zero are given by: 

■ , , m* «£» 1 (t) 

T,(t) = - 7 5 ■■, =“ 

^ Uj(t) + U 2 (t) 

(4.11) 

m* w 2 (t) 

t 2 (t) = - — == r 

/ Wj (t) + w 2 ( t) 


where / T^(t) + T^(t) - w* (4.12) 

Here we conclude that for a symmetrical spacecraft with w 3 = const., and 
T 3 (t) = 0, the control torques required to reduce the transverse 
angular momentum are not explicitly dependent on the type of extension. 
It is doubtful that boom extensions alone could be used to effect time 
optimal control. 

In general, we can say that in. the event of control jet failure 
the booms could certainly be used as a back-up reusable systemfor 
de-tumbling (even if they cannot directly implement time optimal 
recovery) . 


2. Extension of End Masses 

The extension of four end masses, along the symmetry axis of the 
spacecraft, is shown in Fig. 4.1. This scheme is considered as a 
possible means of reducing the transverse angular velocities in a time 
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optimal manner. The end masses are assumed to be equal and they are 

placed very close to the symmetry axis such that d/&j << 1* The 

• * 

control variable are the extension rate £j(t) and T 2 (t). 


The moments of inertia about the principal axes can be written 


as 

Ii - I 2 = I* + 4mU? + l\) = I 

I " l* 

A 3 A 3 


(4.13) 


The equations of motion can be developed, with - w 0 = const, 


as: 


03 1 " 


_ (I* - it) 


03 nO) 


0^2 


^ (4 + 4)(« 0 “2 - “ 1 ) 


- —Jr (t\t\ + T. 2 ^- 2 ) 03} (4.14) 

o> 2 = ~ ~y --- 030031 - (£l + ^2)(wo“l + w 2) 

• * 

- p- (£ 1^1 + ^ 2 ) “2 (4.15) 

Eqs. (4.34) and (4.15) can be rearranged to yield: 

Wj = <i)(i) 2 + g(t)(aij - + ^(0 “i (4.16) 


= “ 0303 ^ + g(t) (<L 2 + WqW^ +h(t) w 2 


(4.17) 


71 



/ T * _ T * \ 

where w - - ^ 

I 


g(t) = - ^ (l\( t) + t 2 ( t) ) 


h(t) = - pr (^(t) ^(t) + l z ( t) c 2 (t) ) 

and we observe, h(t) = g(t) 

Using the transformation, 




COS coot 

-sin wot 


r i 
WJ 

^2 


sin wot 

COS wot 


W2 


Eqs, (4.16) and (4.17) become 


(4.18) 

(4.19) 


(4.20) 

(4.21) 


(4.22) 


sii(t) 

fi 2 (t)j 


0 (ui-u>o) 

-(u-uo) 0 


!2l(t) 

n 2 (t) 


dt 


g(t) 


W 1 


y0 2 (t )j 


(4.23) 


Comparing Eg. (4.23) with tho standard form (Rsf. 8, p. 599, Eg. 


(7.378)) 


x i(t) 


’ 0 

w 


X 

r+ 

j 

+K 

u x (t) 

x (t) 
U 2 J 


-w 

0 


x (t) 
_ 2 


u (t) 

L 2 J 


(4.24) 
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we can write. 


.A (g(t) a At) ) = K u (t) 
dt 1 

_d (g(t) ft 2 (t) ) = K u z (t) 


(4.25) 


For time optimal control 

ui(t) = ±1 


(4.26) 


u 2 (t) = ±1 

Expanding Eq. (4.25), using Eqs. (4.19) and (4.20) with Eq. (4.26), 
the following result is obtained. 


- pF [(£i(t)+£ 2 (t))n 1 (t)+ 2 (£ 1 (t)£ 1 (t)+£ 2 (t)£ 2 (t))n 1 (t)]= ±1 (4.27) 


" [(i.j(t)+^2(t) )fi2(t)+2(£j (t)i;(t)+^. 2 ( t)i 2 (t) )n 2 ( t)]= ±1 (4.28) 

♦ 

From the above equations, we observe that the control variables £j(t) 

» 

and £ 2 (t) are nonl inearl y coupled with the state variables and the 

9 

solution for £j(t) becomes trivial. It can be concluded that the time 

optimal control of this type of system using only boom extension rates 

can not be established by analytic means. However it may be possible 

13 

to consider this problem by using techniques of dynamic programming. 
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FIG. 4.1. EXTENSION OF END MASSES ALONG ’3' AXIS. 
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V. HINGED SYSTEM 


1 . Derivation of Kinetic Energy 

The hinged system to be studied is shown schematically in Fig. 

5.1. The co-ordinate system representation is shown in Fig. 5.2. The 
system consists of a spinning spacecraft with masses attached to mass- 
less booms of constant length l, which in turn are attached to the main 
spacecraft at radius ro. The end masses are released at t = to and 
thereafter swing out from the spin axis. The angles between the booms 
and the spin axis are denoted by a : and a 2 as shown in Fig. 5.2 and 
are initially zero. A special case of this type was considered in 
Ref. 2 (where it was assumed that the transverse angular velocities 
during deployment remained at zero) but here we consider the general 
three dimensional deployment dynamics. 

The development of the kinetic energy of this type of hinged 
system from first principles is considered below: 

The total kinetic energy of the system, in terms of rotational and 
translational energies, can be written as, 

T = T r + Tj- + const, due to (circular) orbit motion (5.1) 

1 2 2 2 

where Tp = ^ (IjWj + I 2 ^2 + I 3 M 3 ) 

T t = 1 M Vn/an + ^ m { V^/cm 

(M = mass of main body) 


(5.2) 
(5. 3) 
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From the definition of center of mass of the system: 


n 

m X V° 

rcm/o - —n < 5 - 4 > 

M + f mj 
i = i 

where point 'o' is the center of coordinate system and the masses are 
assumed to be equal (m-j - m). The velocity of the various components 
relative to the system center of mass may be expressed: 


V roi /cm = V mi /o + V™ (5.5) 

V^/cm = V^/o + V 0 /cm (5.6) 

The components appearing in Eqs. (5.5) and (5.6) can be further repre- 
sented as: 

♦ 

Vnvj/o = ri (5.7) 

V 0 /M = 0 (5.8) 

- m-I.H 

V 0 /cm = - V cm /o = - r cm /o M + ( 5 * 9 ) 


Upon substitution of Eqs. (5.7), (5.8) and (5.9) into Eq. (5.3), the 
translational energy may be expressed as 

T t = £ |V 0 /cm| 2 + 1 l mi |V mi /o| 2 ' 

+ — l mi | Vo/cm| 2 + l mi (Vm-j/o) - (Vo/cm) (5.10) 
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After some algebraic manipulations, we obtain, 

Tt- 7 iCh • v,) - »i • l Vi) (5.11) 

where Vi = ri + wx ri 
M = M + J m-j 

Thus, the total kinetic energy of the system is given by: 

T = \ ( + l 2^2 + l 3*V + f ^ * V 

2 n - n 

- ” (.1 Vi • J Vi) + const. . (5.12) 

2M i" 1 i - i 

As an example, we consider the case from Ref. 2 where m = m/2, 

a 1 - a 2 = a, I 3 = I and w 3 = O. The kinetic energy is then obtained 

as (neglecting orbital motion) 

T = 1 | e 2 + ~ [ l 2 a 2 + e 2 (r 0 + l sin a) 2 ] 


mi i 2 S i n 2 a :2 
2 (M + id) 


(5.13) 


which corresponds identically with Eq. (18) of Ref. 2, which was 
presented without development. 

Next, a more general case of the hinged deployment system con- 
sidered is shown in Fig. 5.3. Here there is no restriction on the 
location of the hinge points. The co-ordinates of the two masses are 
given by 
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(5.14) 


x : = 0 x 2 = 0 

y : = r 0 + l sin a 1 y 2 " -(r 0 + £ sfn 

Z 1 = a * ” ^ C0S a l Z 2 = a * “ ^ C0S 

Here ’a*' is the offset of the hinge point from the '2' axis. Upon 

substituting Eqs. (5.14) into Eg. (5.12), and after algebraic simplifi- 
cations, the resulting equation for kinetic energy is: 

T = \ £ Vf + 1 2 W 2 + Vl 1 

+ ^ [l 2(r 2 +a 2 +£ 2 )+2£ {r 0 (sina 1 +sina 2 )-a^(cosa 1 +cosa 2 ) )} uf 

+ {2a 2 -2a^(cosa 1 +cosa 2 )+£ 2 (cos 2 a 1 +cos 2 a 2 )} 

+ {2r^+2r 0 £(si na 1 +sina 2 )+£ 2 (sin 2 a 1 +sin 2 a 2 ) } 

- {21 {a^(sina 1 -sina 2 )-r 0 (cosa 1 -cosa 2 ) } 

- 2£ 2 (sina 1 cosa 1 - sinot 2 cosa 2 )} w 2 a> 3 

+{2f 2 (a i - cti) + 2 1 {on (ro sinai - a*cosc*i) 

- a 2 (ro sina 2 - a*C 0 Sa 2 )}} 101 + l 2 (aj + a 2 )] 

rc 2 [{ 2 ( 2 a| + l 2 ) + 2 £ 2 cos ( 04 + a 2 )- 4 a*£(cosa 1 +cosa 2 )}iA)f 

” 2(H+2m) 

+' {2a*- l( cosai+cosa 2 )} 2 w 2 + f 2 (sinaj-sina 2 ) 2 
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~2£(sina 1 - sina 2 ) {2a* - ^(cosaj + COSa 2 }} to 2 a) 3 
+Zt{l[a l - a 2 ) + t COS {«! + a 2 )(a! •- a 2 ) 

- 2a* (cOSaj ct j - COSa 2 a 2 )} 

+ l 2 {af + a| - 2a x a 2 COS^ + a 2 )}3 + const. (5.15) 


2. Development of Equations of Fiotion (Neglecting External Torques) 

The equations of motion in the five variables w 1# w 2 , w 3 » a i anc * 

14 

a 2 are developed using the Quasi-Lagrangian formulation for , 
i = 1,2,3 and the general Lagrangian formulation for the variables 
a x , a 2 . The equations of motion for this system can be represented 
by: 


d 

aT 


aT 


aT 

dt 


-U3 3 • 

3o) 2 

3o) 3 

d 

aT 


aT 


aT 

dt 

3co 2 


3o)3 

+ u) 3 

3oji 

J_ 

aT 


aT 

+W 1 

aT 

dt 

3w 3 

— cu 2 

3io j 

3u 2 

d 

aT 

aT 


3F 


dt 

aij 

3a 

— + 
1 



_d 

aT 

- u 

— + 

aF 


dt 

3a 2 

3a 

2 

aa 2 



=* 0 


= 0 


(5.16) 

(5.17) 

(5.18) 

(5.19) 

(5.20) 
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where T = Total kinetic energy of the system 
F = Rayleigh dissipation function. 

Making the approximation: m 2 /M << m or (m/M << 1) and letting F ~ 0 

(for the case of no assumed energy dissipation), the equations of 
motion are obtained as follows: 

Il“r(l2"I 3 )“2 W 3 + m[ 2 ( r 0 +a * + ^ 2 ^ + 21 ^ r 0 ( Sa 1 + Sa 2^“ a ^^ Ca i +CCt 2^ ] ^l 
-m{2a 2 “ 2 a^-f(ccc 1 +ca 2 ) + ^ 2 (c 2 ci 1 +c 2 a 2 )- 2 rg- 2 r 0 -€(sa 1 +sa 2 ) } ^2^3 

+ 2m t { ro(caiai + Ca 2 a 2 ) + a-fc^Sajaj + Sa 2 ct 2 )) 

+ ml {a^ (sa 1 -Sa 2 )-r Q (ca 1 “Ca 2 )- ^ (s2a 1 -s2a 2 )> (« 2 “ ^ 2 ) 

+ ml {l (c^ - a 2 ) + (^i) 2 (^o Ca l + a * Sa l) + ^l( r 0 Sot l" a * Ca l) 

- (a 2 ) 2 ( r o Ca 2 + 3*Sa 2 ) “ a 2 ( r o Sa 2 ” a*Ca 2 )) = 0 (5.21) 

I 2 a) 2 -(l3-I 1 > to 3 u 1 +m{2a 2 -2a^e(ca 1 +Ca 2 )+^ 2 (c 2 a 1 +C 2 a 2 )}w 2 

-mlia-icisa-L- Sa 2 ) “ r 0 (ca a - ca 2 ) - ~ (s2^j - s2a 1 )) a >3 

+ mtCa^Saj - Soc 2 ) - r 0 (cetj “ ca 2 ) - £ (s2aj - S2^. z )} WjWj 

- Rl {l 2 (s 2 a i '+ S 2 a 2 ) - 2 (a 2 + l 2 ) +'2£a + (cai + Ca 2 )) W 3 UJ} 
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+ ra{ 2a*£ (scticti + S<x 2 ^z) ~ -C 2 (s2aia) + s2a 2 a z )} ui 2 

+ m i{ l (cij - a 2 ) -2a jt .(ca 1 a 1 - Ca 2 a 2 ) + l ( c2it 1 a j - c2a 2 a 2 )}u 3 = 0 

(5.22) 

I 3 « 3 - Ol * h )w 1 o ) 2 + m{ 2 r^ + 2 r 0 £ ( sa j + sa 2 ) + Z 2 is ? a 1 +s 2 a 2 )} w 3 
- Sa 2 ) - r 0 (ca : - Ca 2 ) - L (s 2 a 1 - s 2 a 2 )} w 2 
- rn { 2 (r§ +£ 2 ) + 2 .lr 0 (sa } + Sa 2 ) “ Z 2 (c 2 a 3 + c 2 a 2 )} W^W 2 
-me' {a*(sai - $a 2 ) - r 0 (ca x - ca 2 ) - L (s 2 ^ x - s 2 a 2 )} 

-me' { 2 r 0 ($a 1 a 1 - Sa 2 ot 2 ) - l (c 2 ai<*i - c 2 a 2 a 2 ) + l (a 3 - a 2 )} w 2 
+ m { 2 r o e(ca 1 a 1 + Ca 2 a 2 ) + l 2 (s 2 ^ia 1 + s 2 a 2 3 2 )} 0J 3 “ ^ ( 5 . 23 ) 

Z a, + [Z + Sa - cL Ca ) a), 

1 1*11 

2 2 

- ( r 0 Ca l + J s2a l) “3 ~ ( a * Sa l - j s2a l) ^2 

- (r^Ca^ + a*Sotj ) Wj + ( a * Ca 1 + r o Sa i “ ec2a 1 )u 3 &» 2 = 0 (5.24) 
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i 02 - (l + r Q Sa 2 - a*ca 2 ) Uj 

- ( r 0 Ca 2 + f s2a 2 ^ w 3 ~ ( a * Sa 2 " ^- S ^ a 2 ^ U 2 

-(r Ca + a*Sa ) - (a*Ca + r Sa - £c2a 2 )u> tu = 0 (5. 25) 

02 *21 *2 02 32 

where sa-je sina^ and ca-j = cosa-j 
3. Numerical Results 

The five nonlinear equations of motion for the hinged system are 
used to study the torque free motion of the system. The equations have 
been coded for computer simulation and the results are expected in the 
near future. 
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3 


z 


X 2 “ 0 

y 2 = -(r 0 + £ sin a 2 ) 
z 2 = - £ cos «2 


FIG, 5.3, MORE GENERAL CASE OF HINGED DEPLOYMENT SYSTEM. 


xi = 0 

yi = r 0 + l sin cxj 

z\ = a - £ cos otj 

* 
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VI. CONCLUDING COMMENTS 


As a result of the present analysis and numerical results, the 

following conclusions can be made: 

.1. For both' types of telescoping systems, closed form analytical 
solutions for the transverse components of angular velocities as a 
function of time are obtained for the special case where the space- 
craft hub (main part) has a nearly spherical mass distribution and 
where the telescoping system is assumed to originate from the hub 
mass center along one of the transverse axes only. 

2. When the telescoping system is assumed to consist of two 

identical sets of two orthogonally mounted booms in a plane normal to 

the spin axis, the spin axis remains an axis of mass symmetry and, for 

this special case, the analytical solutions are identical to those 

3 

obtained previously. For this special situation it is seen that the 
amplitudes of the transverse components of the angular momentum 
remain constant (but at an accelerating frequency) during deployment. 

3. For the more general case where the hub is not spherical a 
series solution is obtained about t = 0, an ordinary point of the 
time dependent coefficients, in the differential equations of the 
rotational motion of the telescoping system. However, the tadius of 
convergence of such a solution is limited due to the other singular 
points in the coefficients. 

4. The approximate analytical solution for the nearly spherical 
hub and the series solution for the general case are compared with 
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numerical integration results. It is observed that the analytical 
solution corresponds more closely with numerical integration results 
when the extension rate is increased. The series solution can be 
used only in the initial part of the extension where the analytical 
solution also gives essentially the same result. 

5. With fast extension rates and large end masses, the numerical 
study shows that the oscillatory nature of the responses of the trans- 
verse angular velocity components can be reduced rapidly. 

6. As an application for spacecraft rescue and recovery, when 
booms are extended along all the principal axes to detumble a symmet- 
rical spacecraft, exact closed form analytical solutions are obtained 
for all three angular velocities of the spacecraft. 

7. The necessary conditions for asymptotic stability during the 
detumbling sequences can be obtained using Lyapunov's second method. 

The conclusions are that: (1) as time becomes extremely large (and 

boom lengths become infinite) it would be theoretically possible to 
despin a tumbling spacecraft and achieve a zero inertial angular velo- 
city state (of course, such a situation will, in practice, not occur due 
to finite length booms); (2) the final spin about one of the princi- 
pal axes can be achieved by extending all telescoping booms until the 
desired spin rate is reached and then continuing the extension of the 
set of booms along the nominal spin axis until the transverse compo- 
nents of angular velocity reach an acceptably small amplitude. 

8. Numerical examination of other cases for asymmetrical hubs 
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also verifies the practicality of using movable appendages for the 
initial detumbling of randomly spinning spacecraft. 

9. Simple boom extension maneuvers .alone can not be used to 
detumble a randomly spinning spacecraft to achieve a desired final 
state in a time optimal manner. 

10. The constraints on the telescoping system as used for 
detumbling arer (1) the limitations on the extension rate* size of 
end mass masses and the length of booms that are practicable; (2) the 
limitations on the rate of initiable tumble that could be handled by 
the system without compromising its structural integrity. 



VII. FUTURE WORK - FART II 


It is proposed that this effort v/ill be a continuation of the 
research accomplished during the first year (May 1974-May 1975) on the 
dynamics of spin stabilized spacecraft with movable appendages. Part 
I concentrated on the analysis of the motion of a spinning spacecraft 
during the deployment of two types of movable appendages - the tele- 
scoping rod type of varying length during deployment and fixed length 
appendages whose orientation with respect to the main hub can vary. 

In addition the use of these appendages to detumble a spacecraft with 
a random spin to achieve final states of (1) close to zero inertial 
angular rate and (2) a final spin rate about one of the principal axis 
was also considered. In the effort proposed for Part II the following 
will be treated: effect of energy dissipation during deployment; use 

of appendages to detumble spacecraft when the appendages may not be 
deployed along principal body axis of inertia; examination of linear 
optimal control theory as applied to the deployment maneuver by 
selecting different integrand functions in the cost functional; and an 
examination of the effects of first order perturbations such as due to 
solar pressure, gravity-gradient, and small amplitude flexibility of 
the appendages. 

With reference to Table I, the items denoted by an asterisk were 
not treated during the first year, Part I. This Plan of Study has 
been modified slightly in accordance with technical discussions held 
at NASA-Langley. The proposed items for future study, as indicated 
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by asterisks will now be discussed. 

The effect of damping during deployment can be studied by incor- 
porating additional degrees of freedom in the mathematical model of 
both types of appendage systems. For example, a pendulous type of 
nutation damping mechanism on the main hub could be considered and the 
Lagrangian equations of motion for the hinged system modified directly to 
include generalized coordinate(s) associated with the damper motion. 

The EuTerian equations of motion as derived for the telescoping system 
would also require appropriate redevelopment since the center of the 
spacecraft hub would no longer be the instaneous system cente) of mass. 
The effect of energy dissipation during a general deployment maneuver 
could be evaluated using numerical integration techniques. For 
deployment with a small nutation angle - i.e., transverse momentum 
components small when compared with the total momentum, approximate 
analytical approximations such as energy sink will be studied. 

In the area of detumbling (a randomly spinning spacecraft) the 
use of telescoping appendages offset from the principal axes will be 
considered. An attempt will be made to reformulate a modified 
Lyapunov function, either in terms of cross products of inertia using 
the original hub symmetry axes system or in terms of the instantaneous 
principal moments of inertia. Numerical simulation of this more com- 
plicated system will be performed and results compared with those for 
the simpler system. In addition the use of the hinged type system in 
conjunction with a pair of telescoping booms along the 3 axis could 
be examined. 
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The difficulty ini determining a control sequence of extension 
rates for different pairs of telescoping booms which would yield a 
time-optimal recovery of a tumbling spacecraft is seen in Section IV. 
The problem has been that when the equations are written in standard 
state form— e.g. for a case of two sets of booms parallel to the spin 
axis - (where symmetry about this axis is maintained during extension), 
the control function (two different extension rates) is non-linearly 

coupled with the state variables. 

Instead of considering only time optimal contol of a tumbling 
spacecraft, it was suggested at NASA-Langley that linear optimal con- 
trol theory might be applied where now the integrand function in the 
cost functional would contain a quadratic form in the state variables 
plus some function of the control. After appropriate linearization of 
the system an attemptwill be made using the matrix Riccati equation to 

yield solutions for boom extension rates. 

A time optimal control solution for this problem can be obtained 

numerically using the techniques of dynamic programming (gradient 

13 

techniques). This approach was recently employed by Kunciw in 
analyzing the optimal detumbling of the system treated in Ref. 7. A 
dynamic programming solution as applied to the present problem will be 
considered, especially if the application of linear optimal control 
techniques does not yield meaningful physical results. 

As time permits at the end of the study, it is planned to briefly 
examine the effect of such perturbations as gravity-gradient torques. 
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solar pressure, and first order flexibility. It is hoped that this 
effort would establish limits on the rate of tumble that could be 
handled by extendible appendages without compromising their structural 

integrity. 

It is felt that by analyzing the dynamics, control and perturba- 
tions of such types of systems with various types of appendages, a 
valuable insight into the dynamical behavior of more complex systems 
can be obtained. 
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TABLE I - TWO YEAR PLAN OF STUDY 

THE DYNAMICS OF SPIN STABILIZED SPACECRAFT 
WITH MOVABLE APPENDAGES 

CONTENTS 

MOTION DURING DEPLOYMENT 

Spinning spacecraft - small transverse momentum 

1. Hinged Type 

- development of equations of motion 

2. Telescopic Type 

a. End mass moving b. Uniformly; di stributed mass moving 

- Analytical solution for spherical ■■ Hub 

- Series solution for non-spherical Hub 
*3. Effect of Dampers 

USE OF APPENDAGES TO DETUMBLE SPACECRAFT 

1. Telescopic Type 

- derivation of kinetic energy 

a. Achieve zero inertial angular rate 

- Lyapunov Function - Kinetic Energy 

b. Achieve spin about principal axis 

- Lyapunov function - Modified kinetic energy 

*2. Telescoping appendages offset from hub principal axes 
*3. Appendages + "3" axis boom 
OPTIMAL CONTROL 

*1. Aopl i cation of linear optimal 

’control theory using different performance indices 


*2. Use of gradient technique 



EFFECT OF PETURBATIONS 


1. Gravity-gradient 

2. Solar pressure 

3. Flexibility vrith small amplitude 


^Proposed for study in second year (Part II) 
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COMPUTER PROGRAMS 




1. End Hass Moving - Numerical Integration 


PROGRAM FOR G E N . E U L JE R_I_ AN F 0 R_M. U L A T ION 

TELESCOPIC TYPE A. END MASS MOVING 
SUBROUTINE RG SO 1 { T , W , DW ) 

DIMENSION W ( 3 ) , D W ( 3 l_ _ 

REAL II, 12, 13, LT7l 2,L3, 110, 120, I 30, Ml, M2, M3 
COMMON I 1 , 12 , 1 3, 1 10, 1 20, 130, Ml, M2, M3, Cl, C 2, C 3 


Ll-Cl-T 

L2~=C2*T 

L3=C3*T 

DLJ 5 C 1 

DL2 =C 2 
DL 3 = C 3 

A1=_M 1 - L 1_*_L 1 

A2=M2*L2*L2 

A3=M3*L3*L3 

JIM 10 + 2.0*(A2_+A3) 

I 2 = 12 0 + 2 * OM A 3~+~ A 1 ) 

13=130+2. 0*( A1+A2) 

_R 1 = M 1 * 1. 1 -0 L j 

B 2 = M 2"* L2 * DL 2 
B3 = M3 * L 3*0L 3 

_DI_1=A. 0^(02 + 33) 

D I 2 =A . 0£ ( 3 3 + B 1 ) 

DI3=4.0*(B1+B2) 

J^IUr AU 2 r..l ? S 2 3 ( 1 ) ) / 1 1 

D W ( 2 ) = ( l I 3 * “ I 1 ) * W t 3 ) * Vv ( 1 I — 0 1 2 ^ > i ( 2 ) ) / I 2 
DW( 3) = ( ( I 1-12) *WC 1 >*K( 2 ) -0 I 3*w< 3) ) / I 3 

RJJURN 

END 

SUBROUf INE__ RGS02 ( T , W , DW , IHLF,N,P) 
DIMENSION W ( 3 ) , Dw ( 3) , DUMMY! 3) 

REAL 11,12,13 

COMMON 11,12,13 • ___ 

DATA _ DEG/57. 2957795/ 

CALL RG SO 1 { T j K , DUMMY ) 

H 1 = I 1 * J { 1 ) 

. H 2 = 12 *.M 2 ) .. 

H3= I 3*W( 3) 

THETA=ATN2( SORT (Hi*Hl+H2*H2) , h3)*DEG 

TP=T+0. 00005 ___ 

WRI TE (5, 1 ) TP, W, THETA, DW , I MU- 
RETURN 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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1 FORMAT ( 1 X , F9 . 4 , 7F 1 3 • 7 , 110) 

C 

END 

external rg'Yoi7rgso2 

DIMENSION P A R M ( 5 ) , W { 3 ) , DW ( 3 ) , S I IB { 3 ) , WORK { 8 ,3) 

DATA N /_} / 

REAL " 1 1 , 12 , I 3 , I 10 i I 20, 130, M 1 ,M2, M3 

COMMON Il,I2,I3,IlO,l20,l30,Ml,M2,M3 f Cl,C2,C3 


DATA CARDS -- 10 COLUMNS FOR EACH VALUE 

1“ JiLAX,_INn TIAL STEP, TULERENCE 

2- MASSES 

3- INITIAL I*S 

4- C 1 S 

5 - IN I T I A L W * S 

6- TYPICAL SIZES OF W*$ 


RE AD (2, 91) T M A X , S T t P ♦ T 0 L 
RE AD ( 2 i 9 1 ) Ml , M2, M3 
RE AD ( 2^, 91) 110,120,130 

RE AD (2 , 91 ) Cl , C2 ,C3 
R E AD ( 2 , 9 1 } W 
R E AD { 2 , 9 1 ) SIZE 

"WRITE ( 5 ,92 ) t MAX , ST E P , TOL 

WRITE (5, 93) Ml, M2, M3 
WRITE ( 3 ,_9A ) 110,120,130 

WR I TE ( 5"; 95 ) C 1 , C2 , C 3 ' ' 

WRI TE ( 5 , 96 ) W 

_WR I T E ( 5 , 9 7_)__ SJ_Z_E 

WRI TE (‘5798) ' ‘ ’ '* "" 

PARM(1)=0*0 
PARM< 2_) =TMAX 

PARM (3 ) == ST E P ‘ ‘ " ’ ' " - : 

CALL RKSCL !N,S 1 Z E , U W , T OL , P ARM ) 

CALL RKGS l P ARM , W , DW , N , I HLF ,RGSO I , RGS02 , WORK ) 

W R I TE (’5 » 99 j I HL F “ 

CALL EXIT 

9 V F ORMA T ( OF 1 0*0) : 

92 FORMAT! 1 1 T M AX= * , F 8 , 2 , 1 OX , 1 S TE P = * , F 8 . 4 f 1 OX , 'TOL=» ,F8.6) 
9 3 FORMAT ! * OMAS SES * , 3F 1 0 . 6 ) 

94 F 0 RMAT ( ’ U I NIT T * , 3 E 10, 6 1 ~ 

9 5 FORMAT ( * OC 1 , 3F 10. 6 } 

96 FORMAT { ' GW * , 3F 10. 6 ) 

9 7 F CRMA T { ’ 0 S I Z E * » 3 F l 0 . t> ) 

9 8 FORMAT! * 1 • , T 6 , • T • , U 7 , • W 1 • , T 30 , • W2 1 , T A 3 , • W 3 1 , 

_ 2T56, * THETA* , Tb9, *Dwl * , T81 , *DW2 • , 

3T 9 A , • DW3 * , T 1 08 , 1 I HL F 1 , / ) " v 

99 FORMAT ( *OIHLF= * , 13 ) 



ORIGINAL PAGE IS 
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2 . End Mass Moving - Analytical Solution 


C_ ANALY. CAL. W1,W2. A. END MASS MOVING 

REAL I 0 » M 

R E A 0 ( 2 , 5 1 } W 1 0 , W 2 0 

READ ( 2 ,51) DWLO, DV.20 . 

R E A 0 ( 2 , 5 1 ) ToY.M, C 

51 FORMAT ( 5F 16.0) 

fcRlTE(5,52) w 1 0 » W 2 0 

V. R I T E ( 5 » 5 3 ) D V. 1 0 » C W 2 6 
wRITE(5,54) I 0 , M i C 

52 FORMA T { ' lwlO=* ,F15.6, 1_0X, 1 W20 = * , F 1 5 . 6 ) 

•'53 FORMAT {’ 00 W 1 0 = ' , Flt>. 6 , 10X , *UW20^ ' , F15.6) 

5A FORMAT ( ‘ 0 I 0 = * , F8 . A T 1 OX , ' M= • , F 8 . A , 1 OX , * C - * t F8 . A > 

WRITE (5, 5) _ .. .. 

5 FORMAT { » 1 • , T6, 1 f 1 , 1 1 7 , ' Wl ' , T30 , ' W2 * ) 

T = 0 . 0 

STEP = 1 . 0 _ 

15 P = 2 . "6* i 1 * c * c 

A 1 = S £• R T { 2.0-10/P) 

B 1 =SQRT ( 0.5-10/P )_ 

C 1 = T / A 1 

D 1 = A T A - S 1 ( C 1 ) 

E 1 -B 1 *0 1 _ 

F 1 = P/ 10 

G 1 = 1 . 0 + F 1 *T *T 

W1 = U10 + 2.0*DW10*(T-E1) ) / G 1 
A 2 = S C R T ( IO/P) 

B 2 = I 0 / P 

C 2 = T/ A 2 

D2- ATAN ( C2 V * 

E 2 = 02 / A2 
F2=T*T+B2 

G2= T/F 2 

W2 = «20 + 0.5*DW20*R2* (G2*E2 ) 

WRITE (3, 10) T t '.V 1 , W 2 

10 FORMAT ( 1 X i F 9 . A , '2 F 1 3 . 7 } " 

T = H STEP 

IF i T-60.0) 15,15,20 

20 CUM I I MoE 
CALL EXIT 
END 
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3. End Mass Moving - Series Solution 


C SERIES SOLN. CAL. W1,W2. A. END MASS MOVING 

REAL I 10, 120, 1 30, I 1 , 12 , M _ 

RE AD ( 2,21) 110,120, 130,M,C,TF 

READ! 2 , 21) Wl 0 , W20 , DW 1 0 , DW20 , HO 

21 FORMAT (8F 10.0) „ _ 

WR I TE { 5 , 5 1 ) I 10,120 , I 30 , M ,‘C , TF 
HR 1 T E ( 5 • 5 1 ) W 1 0 , W20 , DW 1 0 , UW20 , HO 
51 FORMAT (8F 10.6) 

c ' ' 

P = 2 . 0*M*C*C 

C 

C ' CONSTANTS FDR HI (T j 

AD 1 " I 10* 1 2 0* I 30- I 30* ( 1 20+1 30) 

AE1=I20*(I30*I30*( 1 10+12 0+ I 30 ) +2 . 0*110*1 30 * ( 120+130) )*P 
. A F 1 = i 2 0 * (130 * I 3 6 + I 1 0 * ( I20M30)Y*P*P 
AG 1=120* (11 0+120+3. 0*130) *(P**3.0) 

AH 1 = 1 2 0* ( P * * 4 . 0 ) 

" A1 1 = 2 . 0* I 10* 1*30* I2 6*T2 0*P 
A J 1 = 2 . 0 * ( I 10+130 )*I 20* I 20*P*P 
4K1=2.0*( 120*1 20 )*(P** 3.0) 

A L 1 = ( 1 30-1 1 6 ) ^ < I '30"^ ! '3 0 - 1 2 0^ I 20 ) *HO*HO 

AH 1 = 2. 0* I 30* ( I 30-1 10 ) *P*HO*HO 
AN 1 = ( I .30-1 10)*P*P*H0*H0 _ 

c 

A1=-(AL1 )/{ AD1*2.0*1 .0) 

A2=-( 2.0* 1.0* AEl -2.0* AI 1+ALl )/( ADI *4.0*3, 0) 

A 3 = - A M 1 / ( A 0 1 * A * 0 * 3 . 0 ) 

A4 = - ( A . 0*3.0*AE 1-4. 0*A I 1 + AL l ) / ( AD 1*6. 0*5.0 ) 

A 5 = - ( 2 . 0 * 1 . 0 * A F l - 2 . 0 * A J 1 * A M 1 ) / ( AD 1 *6 . 0*5 . 0 ) 

A6 = -A'U/ ( AD 1*6. 0*5.0) 

A7=-( 6. 0*5. 0*AE 1-6. 0*AI 1+ AL 1 ) / ( AD1*B .0*7 .0 ) 

A 8 = - ( A . 0 *3 . 0* AF 1 - A . 0 * A J 1 + AM ] ) / ( AD 1 * 8 . 0* 7 . 0 ) 

A9 = - ( 2 . 0* 1 . 0* AG 1 -2 . 0* AK 1 + AN l ) / ( ADI *6. 0*7.0) 

A l 6 = - ( 3 . 0*7 . 0* AE 1-6 . 0* A 1 1 + AL 1 > / ( AIH * 10 . 0*9 . 0 ) 

A 11 = ( 6 . 0* 5 . 0* AF 1 -6 . 0 * A J l + AH 1) / ( A D 1 * 1 0 . 0*9 . 0 ) 

A12=(4. 0*3^0* AG 1-4.0* AK1 + AN1)/(AD1*10. 0*9.0) 

A 1 3 = ( 2 .0*1 . 0 * A H 1 )/{ A 01*10. 0*9.0) 

C 

Bl=- (-A fl + ALl ) / ( ADI* 3 . 0*2 . 0 ) 

B2=“( 3. 0*2.0* AC 1-3.0* A I 1 + AL 1 )/( AD 1*5. 0*4.0) 

B3 = -(-AJl+A M l)/( AD 1*5. 0*4.0) 

BA = -(5.0 *4. 0*AF 1-3.0*41 1 + A L 1 )/( AD 1*7. 0*6.0) 
B5=-(3.0*?.0*AF l-3.0*AJl+ AMI )/( AD1*7 .0*6.0) 

B6=— (-AK1+AN1 )/( AD 1*7. 0*6.0) 

B7 = ~( 7; 0*6.0* AE 1-7.0* A I 1 + AL 1 ) / ( AD 1 * 8 . 0* 7 . 0 ) 
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B8=-( 5.0*4.0*AFl-5 .0*A J1+4M1 )/( AD 1*8. 0*7-0) 

B 9 = - { 3 . 0 * 2 . 0 * A G 1 - 3 . 0 * A K 1 + A NJ. )_/ ( [AD 1 * 9 ' . 0 * 8 . 0 ) 

' c ’ 

WR1TE(5,5) 

5 FORMAT ( 'l 1 -, T6 , «JT^ , T_1 7 ’ Wl* t T3 0 t * W 2 ' ) 

T = 0 • 0 * ‘ “ 

S T E P = 1 • 0 

A 1 0 = I 1 0 * W 1 0 ; 

A1 1 = 1 10*DWTO 

15 AA12=A1* ( T**2.0) 

AA1 3= ( A2*Al+ A3 ) * ( T**4 . 0 ) 

‘ " A A 1 4 = ( A 4 * ( i 2 * A 1 + A 3 } +A5*A1+A6 ) * ( T**6. 0 ) 

A151=A7*A4*(A2*A1+A3)+A7*A5*A1+A7*A6 
A152=A8* ( A2-A1+A3) +A9-A1 

AA r 5 = ( A 1 5 n ' A 1 52 ) * ( T **8 . O') 

A161=A10*A7*{ A4*A2*A1+ A3* A4) +A10*A7*A5*A1 
A162=A10*A7*A6+A1 1* ( A4*A2*A1+ 44*A3_+_A5* A 1 + A6 ) 

' A163=A12*(A2*A1 + A3 J+A13 

AA16=(A161+A162+A163)*(T**1C.0) 

A 1 7 = 1 . Q+AA12+AA1 3+ A4 14+ AA1 5+ AA 1 6 

" C 

Bll=Bl*( T**2.0) 

B12=(B2*B1+B3)*(T**4.0) __ 

" B 1 3 - ( B4*TB'2*8 1 + F 3 } + B5*B1 + 36 f* ( T**6 . 0 ) 

B141=B7*B4*(B2*Bl+b3)+B7#B5*Bl+B7*B6 

B142 = B8*(82*_B1+_B3) +09* 8 1 __ 

“ B 1 4 = ( B 1 4 1"+ B142") VCT * * 3 "0T 

B15=1*0+B1 1+B12+313+B14 

C 

‘ ' H1=A10*A17+A1 L*T*B15 
I 1 = I 10+P*T*T 
W1=H1/I 1 


CONSTANTS FOR H2 { T } 

AD2 = I 10*120*130*1 30 

AE2= { 2 . 0* I 1 0* I 20* I 30+120*130*1 30) *P 
AF2=I 20*( I 10 + 2.0*1 30)*P*P 

AG2=I2O*(?**3.0) _ 

A H 2 = 2 . 0* ( 110 + I 30)'**I 2 0* I 30*P 

A I 2 = 2 . 0* ( 1 10+3.0*1 30)* I 20*P*P 
AJ2=4. 0*12 0*(P* *3.0) 

AK2 = ( 130-1 10)* ( I 30- I 20 ) *HO*HO 
AL2 = ( I 30-1 10)*P*H0*H0 

C 1 = - AK 2 / ( AD?* 2 . 0* U 0 ) 

C2 = -( 2.0*1 . 0 * A F 2 + 2 . 0*AH2 + AK2 ) / ( A 02*4. 0*3.0) 
C3=-( AL2 )/( AQ2*4 .0*3.0 ) 

C4 = - ( 4 . 0* 3'. 0*AE2+4 . 0*AH2 + AK2 ) / ( A 02*6. 0*5.0) 
C5=“12. 0*1. 0*AF?+2.0*AI2+AL2)/( A 02*6. 0*5.0) 
C6=-(6.0*5.0*AF2+6.0*AH2+AK2 ) /{ AD2*8.0*7.0) 
' C7 = -(4.0*3.0*AF2 + 4.0*AI 2+A12 ) / { ADZ *6. 0*7.0) 
CB = -{ 2.0*1 ,0*AG2+2 . 0*AJ 2 >/ l A D2 *0 . 0 * 7 . 0 ) 
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c 


0 1 = - ( A H 2 + A K 2 ) / ( A 0 2 * 3. 0*2.0) 

D2=-( 3. 0*2. 0*AE2+3 .0*AH2+AK2 ) / ( AD2*5 .0*4 .0 ) 
03=-{A12+AL2 ) / ( A D 2 * 5 . 0 * 4 . 0 ) 

D 4 = - ( 6 . 0 * 4 . 0 * A E~2 + 5 .0* AH2+ A K2 ) / ( A 02 * 7 . 0 * 6 . 0 ) 


05 = - ( 3.0*2.0*AF2+3.0*AI2+AL2 ) / { AD2 * 7 . 0 * 6 . 0 5 
„ P 6 f “A J 2 / {AD 7^0 0 ) _ 

A20= I 2 0 * W 2 0 

^ A 21 = I 2 0*DW20 

C 1 2 = C 1 * ( T * * 2 . 0 ) 

C 1 3 = { C 2* C 1 + C 3 ) * ( T**4 . 0 ) 

C 1 4 = { C 4 * ( C 2 * C r+ C 3) + C 5 * c 1 j * ( T * * 67 b ) 
C151=C6*C4*(C2*Cl+C3)+C6*C5*Ci 
C 1 5 2 = C 7 * { C 2 * C 1 * C 3 ) + C 8 * C 1 _ 

C 1 5 = ( C 1 5 1 + C 152) * { T * * 8 . 0 ) 

C16=1 . 0 + C12 + C 13 + 0 4 +C 15 
C 

' " 01 1 =b 1 * ( T **2 . o’ ) v 

D12=(D2*D1+D3)*(T**4.0) 
D13=(D4*(D2*D1+D3)+05*D1+D6)*{T**6.0} 

014 = 1 .0 + 01 1 + 012 + 013 

c 

H2=A?0*C 16 + A 2 1*T*D14 

12 = 120 ‘ "" " " ' " • 

W2=H2/I2 

C 

WR I T F ( 5 , 10) T I W 1 , w 2 

10 FORMAT ( lX,F9*4 t 2FI3. 7) 

T = T + $TF_P 

IF(T-TF) 15,15,20"' " 

20 CONTINUE 
CALL EXIT 
END 
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4. Uniformly Distributed Mass Movi ng - NumericalJjite^t^lL^ 


C 

c 



F R C G R A N FCR G E N . _ _ E U L F_R I AN E CR y EL A T I C N 

L M F C R y L Y CIST y ASS MOVING - NL » v E R I C A L 
SUBROUTINE RCSOI(T|ViL’V%) 


c i y e n s i c_n _ _iJ 2 A juQ.h ill 

REAL fl , 12 » 13, L»LF»I10, 120,130 
CCNRCN II, 12,13, 110,120,130, A»CGN,G, 


L , l F 


L=C*T - 

Aa'= (2 . 0/3 . C ) *C E\’* ( ( A-fL ) ** 3 -A*=> 3 ) +2 
PB = 2*C*CEN !!! C S M (A+L }**2-A**2) 


. C-CCN* 


(L 


F-L ) * A* *2 


I 1 = I 1 C + A A 
12 = I 2 C 


I 3 = 1 ?- C + A A 

cn=ee 

C I 2 = C - C 
C I 2 = F? P 

( I 2- 1 3) - h i/'J o L 1 lb. LllA AAA — 

C vT( 2 ) - ( ( 13 -II ) - V. { 3 ) - u ( 1 ) - C i 2 - u ( 2 ) ) / I 2 
CVs { 3) = ( ( I 1 - 1 2 J * w ( 1 ) * K ( 2 ) - C I 3 - a ( 3 ) ) / I 3 


RJJURK 


JENC 


SUBROUTINE _RG S C 2 (_T OK , _l_H L F_, N , P.)_ 

C I N E N 3 I C N W ( 3 ) t C w ( 3 ) » CUN y Y ( 3 ) 

RE AL 11,12,13 

CCyf’C^N 1 1 _ T 1.2 »_I 3 — 

"C A T A 080/57*2957795/ 

C 

- — CALL "RGSCll T ."wTbt i v y Y*) 

F 1 = I 1 *v U ( 1 ) 

F2=l?_*j)(2i : 

' F 3=1 3- .*"( 3 ) 

TF.ETA=AT\2(SCRT(H1*H1 + H2*F2),F3)*CEG 

TP = T4 C . CCCC5 

— j y E ( 5 " i ) T P , V, i T F E T A , CW » 1 H L F 

* RETURN 

— C * 

1 F C R 0 A T (1X»F9.A,7F13.7.I1C) 

" ENC 
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FXTFPNAL RGSC1 , R G $ 0 2 . — -r- : — 

TiTeNS ICTTTrm'T ) Vk'c 2>.DW(3)tSlZE(3>,fcCRK(bf3> 


CAT A N/3/ 

_R E A L 11 , I ? 

CCMMCN I 1 t 1 2 » 1 3 


, 12, I 1C, I_?0, 


C 

C 

C 

c 

c 

c 

c 

c 

c 


13 0 t l , LF 

11C, 120, I 30, A, DEN, C,L, IF 
i fTTrTlTPN S FCR EACP VALLE 


CA1A CAKCS — 1C COLUMNS FCR 

1- TKAX, INITIAL STEP, TCLERENC6 

2- DENSITY, AIR A.C I_L_S Cf SAT E_L LJ LP_L 

“ 3 -“ inTu aTT'S 

A- C , L F 

5- INITIAL W ' S 

V- TYPICAL SIZES CF VS 
TMAX , STEPiTDL 

CEN,_A_ 

Tic, 120, 120 
C ,LF 

sTze”~ 

WRITE (5,92) TMAX, STEF, TCI. 

WRITEJ 5,9 3LJLPTlA 

WR'I TF ( 5 , 94) 110,120,130 

WR1TE(5,95)C,LF 

WRITE ( 3,96) L 

'WRITE (5,9?) SIZE 
W R I T E ( 5 , 9 E ) 

P A RM ( 1 ) =0 * C 


R E A C ( 2 , 9 1 ) 
_R_E A D_( 2 ,91 l 
R E A D ( 2 , 9 I ) 
RE AC ( 2 , 9 1 ) 
R E A C ( 2 , 91) 
'READ (2 ,9 1 ) 


91 

92 
'93 

94 

95 

96 

97 
9 S 


P A R K ( 2 ) = T M A X 
PARM(3)=STEP 

CALL RKSCL(N,S 1 Z E , L : W_ , LP I- f D . A p ^ - 

C A L L R K G S 1 P A R M , W , C W , N » I H L F » R C- S 0 1 , R G S 0 2 , 

WR 1TE (5,99 ) IFLF 

CALL . EX 1 1 

F C R M A T ( 8 F 1 C . C ) , \r 

F C R Y A T ( * 1 TM AX= 1 , f 6 r 2 , 1.P x *. * S I E P “ 1. » F- 8 ‘A* \S: 

FCR MAT ( ♦ C C E N = 1 , F 8 . A , 1 0 X , 1 A - 1 ,F8.9) 

FORMAT ( 'C1NIT I 1 , 3 F 1 C - 6 ) 

FCR v AT ( ’ CC=’ , F h . 2 , 1_CX , _'_LF_= 1 i F 
FORMAT ( 1 OW ' , 2 FiC . 6 ) 

FORMAT { • OS I ZE 1 , 3F1C . 6 ) 

FORMAT ( ' 1 ' i T 6 » ' T ' » 1 Vl.i*.'*} ' 


U C R K ) 


1CX, * TC L 


•THETA* f Tfc9* , Cwl , ,T81, 


, T 3C , 'W2 
• D W 2 % f 


T43 , • V% 2 ' , 


3194 , »CW3' ,T10£ , 1 IFLF * , /) 

99 FCR M A T ( ' C I F L F = * , 1.3 ) 


ENC 
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UN UN' UN 


5. Uniformly Distributed Hass Moving - Analytical Solution 

• c" mct'ly ' c i cry u err y:cnfcpT:~riTr:~rei s~fcving 

PEftL 1C 

R E ft C ( 2 » b H _h.l_C . V_2_C 

reV£(2iY1) U v. 1 C » CY 2 C 
P E ft C { 2 » b 1 ) JCjCtCtTF 

JLi f cr r/ 1 (_L F .i. 6 ?A. J . - 

C CE NS 1 1 Y Cf PCCR 

P R I T E { 5 » 5 2 ) Vn 1 C » Vn 2 C 

_V\ P I T E (Jj » _5 2 ) 

Yr I tT( 5 1 CC ) I c , c , c » f F 
2 F C R Kft 1 ( 1 U1C= ' , F 1 5 . ft 1 1 0 X , ' ;% ? C = ' » F 1 5 . ft ) 

a F r R y i T ( ‘CC'a1C_=1 u F 15_.6 , 1CX, , .nw2C=_ , .f F .15.6 )._ 

'i'lc^MiY'Cic’Y' , F 8 . A » 1 C X t 1 C = 1 , F £ . A f 1 C X i 1 C = ' »FE.h 
Vn R I T E ( 5 , 5 ) 

5 F C R y A T ( 1.6 ^ 1 T * > T 1 7 i 1 V. 1 _* X .3 j? t.J b_?— 1 1 J 

7 = c7c~ 

5 1 EP= 1 .C 

E = ( { 1 . 5^ I C / C J.J *.Jr -- /-£ — — .. 

U? 1 T t ( 5 . ft 1 1 E 
6 1 F C R R A T ( 1 X , F <= . A ) 

15 n = v. _i c / ( i ._c t j 3..!.? J- A? * °_L 

E xTtf ft’ TERR S in IT F PI ft R E LEFT 
El = ( (0*'- 3.0*1 )/{0**2.CMT** 3.0 

F = { { P + T ) $ * 2 . C ) /_ ( BJ 7 ; 2_ ._C - E * T + T * *_2_. G_ ) 

Cf= ie/3 * C Y* ft L C G ( F ) 

C=22.C/A2.C 

c 1 =_i 2 . c - t - e ) / ;u . 7 ^. 2 i 
c 2 = ft 1 ft ;< ( c i")Y c 
C2=(2.C*B)/1.732 

CA = C 2 -C 2 

E l" = C Vn 2 C / 3 . C 
V 2 = W2C +E 1* ( B 1 +C 1 +CA) 

V P 1 1 1 l *5 1 1 C ) J r V. 1 , > 2 _____ 
i c F C R v ft 1 ( 1 X i'F 9 . A » 2 F 1 3 . / ) 

1 = n 5 t ep 

IF t T- I F ) lb , 1 5 , 2_C 

2C CCM 1 NLE 
C ft L L EXIT 

jiU __ 
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6. Detumb1ing--To Achieve Final Spin Along On e of the Principa} 
Axes (End Hass Moving - Numerical Integration) 


A 

; C 
; C 
; C 
; C 
; C 
; C 
; C 
; C 
; C 
C 


C 


C 


EXTERNAL RGS81,RG502 

DIMENSION PARM(5),U(3) . DU(3) , SIZE (3) ,IJ0RK(8.^) 
REAL 11, J2, 13, 110, 120 - 130, Ml , M2, M3 ^ 

COMMON 11,12, 13. 1 10. 1 20. I30.M1.M2,M3.C1,C2.C3 


DATA CARDS — 10 COLUMNS FOR EACH VALUE 

i_ TI-1AX, INITIAL STEP, TOLERENCE 

2- HASSES 

3- INITIAL I'S 

4- C'S 

5- INITIAL ITS 

G- TYPICAL SIZES OF U'S 


CALL IN0UTC2. 5) 

N- 3 

TYPE ' RKGS JOB' 

FnEAD (2. 9 I ) TMAX, STEP , TO L 
READ (2, 9 1 1 Ml, M2, M3 
REAIK2.91) 110.120*130 
READ (2,91) C1,C2,C3 
READ (2, 9 H U 
READ (2, 9 1 ) SIZE 
UR I TE < 5 > 92 ) TMAX, STEP , TOL 
WRITE (5, 93) Ml, M2. M3 
UR I TE (5, 94) 110,120.130 
WRITE (5. 95) C1.C2.C3 
WRITE (5,36) U 
WRITE (5. 97) SIZE 
PARM(l) =0.O 
PARMC2) =T!1AX 
PARMC3) =STEP 

CALL RKSCL ( N , S I ZE . DU, TOL , PARM) 


WRITE (5, 98) Mnot „ 

CALL RKGS C PARM, W, DU, N, IHLF, RGS01 ,Rlso0«-,U0RK) 

UR 1 TEC 5, 93) I HLF 


CALL EXIT 


91 FORMAT ( SF19.0 ) 

92 FORMAT ( ' 1TMSX=' ,FS.2, 10X, 
33 FORMAT ( ’ OMASSES' , 3F 10 . b) 

94 FORMAT ( ' O IN I T l',3FlG.b) 

95 FC'RMATc'VC ',3F10.5) 

FORMAT ( * CUJ ' , 3F 1 0 . b) 

9? FORMAT (' OS tZE ',3FlQ.b) 

98 FORMAT ( ' l',Tb,'T',T17.'Ul* 
2T55, 'THETA' ,Tb9. 'DW1' . T81. 
3T94, ' DW3 ' , T10S. ' IHLF',/) 

99 FORMAT ( ' 0 1HLF = ' , 13) 


'STEP 55 '.. F8.4, 10X,'TOL« 


, T30, 'W2' , T43, ' U3' , 
' DW2' , 


,F8.6) 


END 
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C PROGRAM FOR GEN. EULER I AH FORMULATION 
C TELESCOPIC TYPE 0. END MASS MOVING DETUMBL ING 
SUBROUTINE RGS0 1 ( T, U. DU) 

DIMENSION U(3) , DUC3) 

RE PL II, 12, 13. LI .L2.L3, I 10. 120. 130. Ml, M2, MS 
COMMON 11,12, IS, I j 0 > 120, 130, Ml , M2, M3, C 1 , C2. C3 
C 

c 

IF (T.GT. 2 .5) GO TO 20 
L1=C1*T 
L2=C2*T 
DL I -C 1 
DL'2=C2 
GO TO 30 
C 

20 Ll=Cl*2.S 
L2=C2*2.5 
DL 1=0.0 
, DL2 =0 . 0 
30 CONTINUE 
L3=C3*T 
DL3*C3 
A 1 =M 1 *L 1 I 
A2=M2*L2*L2 
A3=M3*L3*L3 
11*M0+2.8*(A2+A3) 

12* 120+2. 0*(A3+fil) 

13 = 1 30+2.. 0H- ( A 1 +A2 ) 

B1*M1*L1*DL1 
B2 = M2*L 2 -DDL 2 
B3«M3*L3*DL3 
Dll =4. G*<B2+G3) 

DI2“4.0tt(B3+Bl) 

D 1 3=4 . Ci'K(B 1 +B2) 

PlJCl) = ( C 12- I3)*U(2)*U(3)-D1 1-U.K 1) )/ \ 1 
DU ( 2 ) = ( ( 1 3- 1 1 ) *U ( 3 ) *IJ C 1 ) -D 1 2*U ( 2 ? ) / 1 2 
DUC3) = C ( 1 1- 12) *UC I ) *U<2) -D 1 3 R.K3) ) /I3 
RETURN 
END 

PARAMETER PEG=57. 2S57795 
SUBROUTINE RGS02CT, U, DU, IHLF, N, P) 

' I) 1 ME NS I ON U C 3 ) , PUC 3 ) , DUMM Y ( 3 ) 

REAL 11,12,13 
COMMON 11,12,13 
C 

c 

CALL RGSO 1 ( T , U, BUMi 1Y) 

H 1 - 1 1 *U ( 1 ) 

H2= 12*LK2> 

H3= 13-1LU3) 

THETA =ATAH2 (SORTCH I *H \ +H2*H2> .H3) *PEG 
TP =T» 0.00005 

UR ITE (5,1) TP, U, THETA, DU, IHLF 
RETURN 

1 FORMAT (IX, F9.4,7F13.7, 110) 

END 
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non 


7. Subroutine RKSCL 


SUBROUTINE RKSCL (NU, UMAC, DU, TOL, P) 
REAL UMAGC1K DU(l), PC4) 


H - HU 
UHORil = 0.0 
DO 1 1*1. H 

UMGRM = UNORM + l.O^UnfiG(I) 
J CONTINUE 

UNORi-1 = J . B^UHORM 
DO 2 1=1. N 

DU Cl; = UNORM/UMAGd) 

2 ' CONTINUE 

PC4) = H*UiiORM*TOL/15.0 
RETURN 
C 
C 

END 
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8. Subroutine RKGS 


SUBF’f'UT I tit RKGS ( PRMT , Y,DEEY»NpI M : I HLF ^ t- C > ' ‘ p ' - ■ 

DIMENSION Ya^DSRYUj.HUAvd, u.H^.bC^^ 

DO 1 I = 1 MID IN 

1 AUX(8, I) s .GbbSbSo<'vi'hHT(U 

X=PRMT(1) 

XEHD=PRHT<2) 

H=PRnt (3) 

PRMTC5) =0. 

CALL FCT(X> Y,DtRY) 

ERROR iESi 

IF (H*(XEND-X) ) 38. 3K.. 4_ K ^, nTi 

PREPARATIONS FOR RUN'-t-KJ i ih IfcinOD 

2 A(i)=.5 

A C2) = . 2928932 
A (3) = 1 . 7u? 10i' 

A (4) = . 1 bbbbb f - 

B(2) = i . 

B (3) * 1 - 
BC4) -2 . 

C (1 ) = - 5 
C<2) = .2928332 
C (3) =1-70? IQ? 

PREPARATIONS OF FIRST RUNGE-KU i ift S i EP 
DO 3 I = 1 > NT IH 
AUXC 1 * I) =Y( I ) 

AUX(2, n=I l ERY(l) 

AUXC3, I)-0. 

3 AUXCh, I)=0. 


RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 

RKGS 


1 

MB 1 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 
. 13 

14 

15 
IS 
17 
IS 
15 
20 
21 
22 



IREC'B 

RKGS 


H=H+H 

RKGS 


3HLF=- 1 

RKGS 


]STEP=0 

RKGS 


1END=0 „ , n 

RKGS 


START OF ft RUNGE-KUuA SitP 

RKGS 

4 

1 F ( ( X+H-XtND) *H ) 7 , 6 , 5 

RKGS 

5 

H*XEND-X 

RKGS 

S 

RECORDING OF INITIAL jWLUES OF THIS STEP 

RKGS 

RKGS 

7 

CAL L OU i P C X, Y, PERT , I R tC, Niu l D i- xh i > 

RKGS 


IF (PRHT C5) )40, S* 40 

RKGS 

8 

1 TEST=0 

RKGS 

9 

ISTEP* I ST tP+1 i nno 

RKGS 


START OF INNERMOST RUNGE-kUTTA LOO. 

RKGS 


24 

25 

26 
2 ( 
28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
^>S 

39 

40 

41 

42 

43 

44 

45 

46 


>1 
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. 10 AJ=A(J) 

CJ=C(J) ? 

DO I 1 I' UNDIM ■ 

RI*H*DERY(I> ' . 

R2 = AJ*<F:1-BjVAUX(6» I ) ) 

Y(I) =Y( D+R2 
R2=R2+R2- L R2 

11 AUX ( 6 , i)=AUX(S> I)+R2-CJ*R1 
IF(J-4)12, 15, 15 

12 J = 

1 h C J ** 3 ) 1 » U' lb 

13 X=X+.5*H 

14 CALL FCTCX, Y, DERY) 

GOTO 10 

C END OF INNERMOST RUNGE-KUTTA LOOP 

C TEST OF ACCURACY 

15 IF C I TEST) I b y IS - 20 _ 

C IN CASE I TEST =9 THERE IS NO POSSIBILITY FOR TESiING Ot- hCCUkhCY 

lb DO 17 I=!.ND:M 

17 AUXC4, i)-Y(I) 

I TEST - 1 . 

ISTEP= ISTEF+IS i tr-2 

18 1HLF=IHLF+1 
X=X~H 

H = . 5*Fi 

DO 19 I - I . ND IM 
Y( 1 ) =AUXC 1,1) 

DE R y C I ) = hLU C 2 , I ) 

19 AUXCb, I) =AUXC3, D 
GOTO 9 

C IN CASE ITEST= 1 ilSiING Oh ACCURACY IS POSSIBLE 

20 I MOD » I STEP/2 

IF C I STEF— I MOD- 1 HOD) 2 1,23,21 

21 CALL FCTCX,Y,DERY) 

DU 22 I»1,HDIM 
AUXC5, I) =Y( I ) 

22 AUXC7, I) =DERYC I ) 

GO 1 0 9 

C CGi'PUTATIGN Or TEST VALUE DELT 

23 DELT=0. 

- DO 24 1=1, ND iH 

24 DELT=DtLT-i-AUXC3, I >*hE!S (AUXC4, 1) -YC1 ) ) 

] F CD E L T- P R KT ( 4 ) ) 2S , 2 8 . 25 

C ERROR IS TOO GREAT 

25 IF C IHLF- 10)26 • 36, 36 

26 DO 27 1 = 1, ND 1 71 

27 AUXC4, i) =huXo, I ) 

I STEP = ISTEP+i 3 1 EP-4 

X=X-H 

IEHD-D 

GOTO IS 

C RESULT VALUES ARE GOOD 

28 CALL FC7CX, Y, DEFY) 


RKGS 

47 

RKGS 

43 

RKGS 

43 

RKGS 

50 

RKGS 

51 

RKGS 

52 

RKGS 

53 

RKGS 

54 

RKGS 

5d 

RKGS 

55 

RKGS 

57 

RKGS 

5 b 

RKGS 

53 

RKGS 

S3 

RKGS 

£1 

RKGS 

62 

RKGS 

6 b 

RKGS 

64 

RKGS 

b“i 

RKGS 

65 

RKGS 

Sr 

RKGS 

6S 

RKGS 

63 

RKGS 

70 

RKGS 

71 

RKGS 

72 

RKGS 

73 

RKGS 

74 

RKGS 

75 

RKGS 

76 

RKGS 

77 

RKGS 

78 

RKGS 

7? 

RKGS 

SO 

RKGS 

81 

RKGS 

G“ 

^ w. 

RKGS 

83 

RKGS 

84 

RKGS 

85 

RKGS 

86 

RKGS 

87 

RKGS 

83 

RKGS 

89 
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l»0 29 IH.NDIM 
AUX<1, I ) -Y( I ) 
AUXC2, I ) =DERY( I ) 
AUX ( 3 . I) = A'jX Lb, I) 
Y< I ) =AUX(5- I) 


29 DERY( I > -hUX(7/ I > 

CALL OUTPCX-H, Y,DERY, IHLF,NDIM, PkMi ) 


IF (PRHT(5) ) 43,30,43 
30 DO 3i I = i -HD In 
Y( I) =hUX< 1,1) 


31 DERYC I) =hUX(2, I) 


IkEC=IHLF 
IF( I END) 32, 32, 39 
INCREMENT GETS DOUBLED 
32 IHLF-IHLF-1 
JSTtP= ISTEP^‘2 


H=H+H 

IF( iHLF)4,33,33 

33 I MOB- I STER/2 

IF< iSTEP- It'iQD-IM0D)4,34, 4 

34 IF (DELT- . 82*FRMi (4))b5,b5,4 

35 1HLF= IHLF- 1 
ISTEP=ISTEP/2 
H=H+H 

GOTO 4 

RETURNS TO CALLING PROGRAM 


30 1HLF= 1 i 

CALL FCTCX, Y,DERY) 

GOTO 39 
. 37 IiiLF = 12 
GOTO 39 

38 1HLF=13 

39 CALL OUTP(X,Y,DERY, IHLF,ND1M,PRHT) 

40 RETURN 
END 
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